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Editorial 


M any Eureka editors have apologised for a delayed pu- 
blication in this space, although most could not blame 
a global emergency One who could have was P. Bro- 
detsky in 1944 - almost two years after the previous issue had 
promised a biannual publication rhythm - but he had the grace to 
not so much as allude to the War, instead blaming the “paucity of 
contributions”. Our publication dates are since so reliably missed 
that we have taken to dehning Eureka as a quasquiannual publi- 
cation wherever we keep a pretence of professionalism, and this 
is almost perfectly accurate: It was published 65 times between 
1939 and 2020, or about once every 1.246 years. I am thus forced 
to admit that the delay of the present issue, while not unusual as 
such, is decidedly above average. 

Eor all its dithculties, we are glad that Eureka is still active 
after so many years, not a matter of course in the volatile envi- 
ronment of student societies. It has certainly seen changes: When 
the hrst issue was published at the very modest price of 6d, it was 
to give students “a medium for stating their views, for discussing 
their present training and their hiture prospects, [and] for pu- 
blishing their less orthodox or less mature researches.” 

The publishers made good on that promise: Eureka No. 2 
had a questionnaire regarding lecturing techniques and the Tri- 
pos syllabus. Besides alterations to the Part II syllabus, results in- 
cluded compelling lecturers to print copies of their lecture notes, 
certainly something that would appeal to students today. The fol- 
lowing issue then called for signatures on a resolution “against 
the calling-up [for warj of students in the middle of their courses, 
because we feel that this will lower educational standards [...].” A 
joint protest to the War Othce was already in planning. 

It was hirthermore reported that “owing to the outbreak of 
the war it has not been possible to make much progress with 
[getting submissions from] foreign students.” Indeed, Eureka had 
always featured articles not only by students and academics from 
Cambridge, but also from other British universities and the wider 
world. The exact balance to strike was the topic of some debate 
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over the years: Issue No. 33 opened with the remark that Eureka had “lost almost 
all its bias towards Cambridge ’ whereas once it was “written by Cambridge men 
and providing a forum for discussion of Cambridges own peculiar problems. 
[...] This could, of course, change overnight with the inauguration of a traditio- 
nalist as editor and [the editor hasj no doubt that this breath of fresh air owes a 
lot to a succession of Churchill editors in recent years.” Eight years later, in refe- 
rence to those Churchill editors, “[thej present editors [...], both Johnians, have 
sought to revive the old values and once again instil the spirit of Cambridge into 
Eureka: all the articles this year have been written by people who have studied at 
Cambridge.” Ear from such extremes, it has since been our good practice to have 
a healthy mix between the two, besides others with an article by Archimedeans 
member Jessica Wang from the University of Durham and a guest article by Josef 
Greilhuber from the University of Yienna. We see this as playing a part in esta- 
blishing links between mathematics students across universities, which is one of 
the stated purposes of the Archimedeans. 

Issue 37 was the hrst to indirectly aid the Archimedeahs other goals, by 
turning a proht. This was due to the introduction of computer topography to the 
editing process, an advantage since all but annulled by subsequent advancements 
such as colour print and glossy paper. Nonetheless, Eureka still contributes to 
the Archimedeahs budget, much to the joy of all our garden party guests. In fact, 
you, too, can make a proht from Eureka: Of the problems posed by Paul Erdos in 
Issue 51, many are yet unsolved, and the cash prizes are still there to be claimed. 
But where on Earth might you turn up a copy of Eureka 51? And why indeed am 
I going on about this magazine’s publication history? 

It’s because all this and much more can now be discovered in our new ar- 
chive, available at 


archim.org. uk/eureka/archive. 

Whether you want to read G. H. Hardy’s famous essay on Mathematics in 
war-time in the original, or the various humorous and serious works published 
over many years by a group of students under the pseudonym Blanche Des- 
cartes, or get a glimpse of Tripos history, it is all there. Besides, while most of 
our articles range from informative to entertaining - with occasional lyrical es- 
capades bordering on the absurd -, a number of original research articles have 
also appeared in Eureka over the years, some of them frequently cited, and we 
felt that these should be made freely available. I am highly gratehil to Michael 
Grayling, Editor of Eureka 65, and Joseph Myers, our archivist, for their work 
with scanning the old issues. 

Eor making the present issue, I thank the contributing authors, and all our 
assistant editors. I am gratehil in particular to Hamza, who brought extraordi- 
nary dedication to the task as an external society member, and to Gesa, who was 
available for some last-minute editing, preventing a hirther delay. With that, I 
hope that you enjoy reading this magazine, and optimistically promise that the 
next issue will hnd its way to you soon. 


Valentin Hilbner 
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Year in Review 

James Bolton, President of 2019-20 


T he past year has been a great year for the 
Archimedeans. We welcomed over 200 new 
members to the society at the freshers fair 
and introductory talk. It is deeply regrettable that 
such a successhil year was cut short by unfortunate 
global events; Tm sure many of our members are 
as saddened as I am at having to forego our annual 
garden party. However, regardless of world goings- 
on, we must - and shall - continue to provide 
mathematical talks and events for our members. 

Our weekly Priday talks throughout the year 
have covered topics from fluid dynamics to sta- 
tistics and black holes to algebraic geometry. 
Highlights included a highly entertaining talk on 
dimensional analysis from Prof. Huppert, a crash 
course in fluid mechanics from Prof. Lauga as well 
as a collaborative talk with the Trinity Mathemati- 
cal Society by Prof. Birkar which was (as might be 
expected) highly popular. 

In addition to our academic events we also 
held a variety of social events. Our ever-popular 
board games and pizza nights proved to be a suc- 


cess drawing many to the CMS throughout the 
year. A collaborative formal with Scisoc and Biosoc 
drew 50 guests. We also introduced new maths and 
art sessions for the first time. The annual dinner 
proved to be the most signihcant social event of 
the year, held at the Royal Cambridge Hotel, it was 
enjoyed by all; my thanks go to Prof. Huppert for 
being our after-dinner speaker. 

Despite not being in Cambridge in person 
in Easter term the Archimedeans continued with 
a vibrant schedule of events. To fill the gap many 
students stepped in and provided talks conducted 
virtually, as well as a virtual pub quiz and integra- 
tion bee. 

My thanks go out to all of our speakers, both 
academic and student, as well as the committee 
and the Eureka team. I would also like to thank our 
members for making the Archimedeans the truly 
great society that it is. I wish all the best to the new 
committee and have no doubts that they will conti- 
nue the good work of the Archimedeans in the co- 
ming year. 
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The Committee of 2019-2020 


President 

James Bolton (Pembroke) 

Vice President 

Jiazheng Zhu (Pembroke) 

Treasurer 

FongTszLo (Magdalene) 

Events Managers 

Beth Holmes (Murray Edwards) 
Shavindra Jayasekera (Trinity) 


Secretaries 

AndyZhao (Trinity) 

Matthew Wearden (St Catharine's) 

Publicity Officer 

Parth Shimpi (GonvilleandCaius) 

Sponsorship Officer 

Michael Norman (Pembroke) 

Webmaster 

Yalentin Htibner (Pembroke) 







◄ 

Those members 
of the then newly 
elected 2020-21 
committe who 
bothered to stay 
for the group 
photo 


► 

Our 2019 garden 
party. Note that 
physical closeness 
was still considered 
acceptable at that 
time. 




◄ 

The Archimedeans' 
committe elections 
have long enjoyed a 
reputation as the most 
rigorous exercise of 
democracy in the world 
of student mathe- 
matics. 


















Student society admits ► 
actual supplies are Anite after 
promising "unlimited ice cream" 
in advertisement slammed as 
"misleading"on Mathmo Memes: 
calls for event managers to resign 


► 

Even a broken leg 
could not stop Segev 
from attending our 
garden party, the 
highlight of his aca- 
demic year. 


















A We were honoured to have 
Prof. Herbert Huppert give the 
speech at the 2020 Archime- 
deans' Dinner. 


Our annual dinner is attended 
T by students from all disciplines. 












Three generations of Presidents: each of them 
enriched the society with their unigue personality. 
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Tetrahedral 
Symmetry and 
High Temperature 
Superconductmty 

Piper Powler-Wright 


The origin of high-temperature superconductmty 
remains a major open problem in condensed matter 
physics. A simple lattice model believed to capture 
the essential physics of the phenomenon is the Hub- 
bard model. In this article we aim to solve the Hub- 
bard model for various numbers of spins on a single 
tetrahedron. Central to our strategy will be the use of 
symmetries. As such, we hrstly review the symmetry 
group relevant to the problem at hand - that of the 
tetrahedron. 


tetrahedron about the three coordinate axes shown in 
Pigure 1. Labelling the vertices 1 through 4 and per- 
forming these rotations, we obtain a representation 
{I^Rx^Ry^Rz} of the group where / is the identity 
matrix and 


Rx 


'0 0 1 0 ' 
0 0 0 1 
10 0 0 
0 10 0 


Tetrahedral Symmetry 

The symmetry group of the regular tetrahedron is 
the collection of permutations of four objects. In 
fact, we cannot make use of the entire group, since 
in general its elements do not commute, a property 
necessary for the simultaneous diagonalisation our 
method will rely on. Instead we turn to the largest 
abelian subgroups of 84 ^. These are (T4, the cyclic 
group of order four, and T^, the Klein-4 group. We 
choose to use ¥ 4 , but C 4 would be equally effective in 
reaching a solution. 

The elements of ¥4 correspond to tt rotations of the 


R 


y 


'0 1 0 0 ' 
10 0 0 
0 0 0 1 
0 0 10 


Rz 


'0 0 0 1 ' 
0 0 10 
0 10 0 
10 0 0 


We will have need of the eigensystem of these ma- 
trices. This is particularly simple as /?^ = / (a G 
{x^y^z}) restricts the possible eigenvalues to ±1. 
Solving the 4-by-4 problem one obtains 
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T 


~-l 

1 


-1 

1 


1 

_ 1 _ 


_ 1 _ 


Rx = +1 Rx = ~l 

Ry = +1 Ry = +1 

Rz — RxRy = +1 Rz — ~1 


■ 1 ■ 


~-l 

-1 


1 

1 


1 

-1 


-1 


Rx — +1 Rx — ~1 

Ry = -l Ry = -1 

Rz = ~1 Rz = +1 


These vectors, as ordered, will be relerred to as the 
Ra basis. 


The Hubbard Model 

The Hubbard Hamiltonian describes the behaviour 
of spins (random variables with two possible values, 
a =^0^4.). It is written by the physicists as 


H=-t + (*) 

(jdO.cT j 

where j labels the sites of the lattice and t,U are pos- 
itive parameters obeying U t. and Cja- are cre- 
ation and annihilation operators for a spin-1/2 par- 
ticle at j satisfying the anti-commutation relations 

{Cjcr’ ^j'cr'} — ^jj'^aa'y 2.nd Tlja \= 

Portunately, one does not need to be familiar with 
these notations in order to work with the model: the 
terms in (*) are simply a set of rules for the opera- 
tion of H on any vector or ket \v) representing the 
physical state of the system. These rules are readily 
understood by example. Take a single up spin on the 
tetrahedron. Four possible states are immediately ev- 
ident: 

2 4 

.P, 

|1> |2) |3) |4) 



These states constitute an orthonormal basis for the 
single spin system. 



A Pigure 1: A tetrahedral arrangement 
of lattice sites and a two dimensional 
schematic (top right) we use to represent 
this geometry. Tetrahedral symmetry: The 
elements of describe tt rotations about 
three axes passing through the centre of 
each pair of opposite faces of the cube 
shown. For example, R^ is the product of 
the two-cycles (13) and (24), expiaining its 
matrix elements on the previous page. 


The lirst term of H, known as the hopping term, de- 
scribes the transfer of the spin from one site j to a 
neighbouring site j'. For instance, 

H\l) = -t\2)-t\3)-t\4), 

since from the bottom left site the spin can hop to 
any of the other three sites. The sum J2{j is taken 
over all nearest neighbours once; i.e., the 6 edges of 
our schematic. 

The second term of H is only relevant in the presence 
of multiple spins: rijj- counts the number of ups at j 
and Uj^ the number of downs. Consider the states 
from a system of two spins: 

| 1 > | 2 ) 

One has 

j 

and 

U Uj^Uj^ |2) = 0. 
j 
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In other words, this is a diagonal term giving -\-U in 
the case of double occupation and zero otherwise. It 
represents an on-site repulsion between, for example, 
negatively charged electrons. In general, any site may 
be occupied by an up, a down, both or neither. Note 
that two spins of same orientation cannot share a site, 
a fact known as the Pauli exclusion principle. 

Using these simple rules we can construct a matrix 
representation of H. Returning to the single spin sys- 
tem. 


ll) 

|2> 

|3) 

|4) 

0 

— t 

— t 

-t 

-t 

0 

— t 

-t 

—t 

—t 

0 

—t 

-t 

—t 

—t 

0 


‘SoIving’ the Hubbard Hamiltonian amounts to soIv- 
ing the eigenvalue problem associated with this ma- 
trix. Of particular interest is the eigenvector with 
the smallest eigenvalue, since this describes the min- 
imum energy or ground state of the system in ques- 
tion. 

In the present case we notice that 


Opposite Spins 

Each spin is free to occupy any of the four sites inde- 
pendently of the other, giving 4^ = 16 basis states, 
four of which invoIve double occupation. Although 
H is no longer of the same dimension as the Ra, we 
crucially use the symmetries to guide the labelling of 
our starting basis in the foIIowing way: Starting from 
an arbitrary state 11), apply each rotation to generate 
the next three states 12), 13), 14). Pick an unused state 
|5) anddo the sameto obtain |6), |7), |8). Repeatthis 
procedure until all 16 states have been accounted for 
(Eigure 2). 



H — —t(Rx + Ry + Rz)- 


Therefore H is diagonal in the R^ basis: 


-3t 0 0 0' 

0 t 0 0 

0 0 t 0 

0 0 0 f 


We conclude that the ground state is non-degenerate 
and of energy —3t. It is characterised by the vector 
[1,1,1,1]^, an equally-weighted superposition of all 
four states, as one would expect on the grounds of 
symmetry. Note that exactly the same result would 
have been found for a single fixed down spin and the 
corresponding set of states- another symmetry. 

This illustrates the use of symmetries in the most ba- 
sic case. In the main part of this article we soIve H for 
two opposite spins, in the limit U ^ t prescribed by 
the model. Systems of 3 or more spins can be handled 
in the same way. 


A Pigure 2: Opposing spins. Successive 
states in each row are obtained by rotat- 
ing the first state in that row tt radians 
around the x, y and 2 : axes in Pigure 1. 
Technical point: The operators cj,| 
anticomnnute, so one should be mind- 
ful of the order in which they occur. We 
choose to always write the up spin first (e.g. 

|1) := 4i4t ^ 

That way, there are never any sign 
changes (+t terms) due to reorder- 
ing after a hop between states 
(Exercise: anyone who has seen second 
guantisation may like to think why this is). 


Writing down H is then a matter of determining 
which states are connected by a single hop of either 
spin; each such connection conveys a — t entry to the 
matrix in accordance with the above rules. Eor in- 
stance, one clearly has (5|i7|l) = —tbut (2|i7|l) = 
0. Since the rotations move two spins, the states in 
any one row are always more than one hop away from 
each other. Remembering that states |13) - |16) also 
come with a U self term, we find 


























16 


u 

u 

u 

u 

1 

—l 

u 

u 

—l 

1 

—t 

—t 

u 

u 

1 

—t 

u 

—t 

u 

0 

0 

0 

0 

1 

1 

0 

—t 

—t 

0 

1 

1 

-t 

—t 

0 

0 

1 

1 

0 

—t 

0 

-t 

0 

0 

0 

0 

1 

0 

—t 

—t 

0 

1 

0 

0 

—t 

—t 

1 

—t 

0 

—t 

0 

0 

0 

0 

0 

1 

—t 

0 

0 

-t 

1 

0 

0 

-t 

-t 

1 

0 

—t 

0 

-t 

—t 

0 

0 

—t 

T 

1 

o“ 

0 

0 

~0 

T 

1 

—t 

0 

—t 

~0 

T 

1 

—t 

—t 

0 

“o 

0 

—t 

—t 

0 

1 

0 

0 

0 

0 

1 

0 

—t 

0 

-t 

1 

-t 

—t 

0 

0 

0 

—t 

—t 

0 

1 

0 

0 

0 

0 

1 

—t 

0 

—t 

0 

1 

0 

0 

—t 

—t 

—t 

0 

0 

—t 

1 

1 

0 

0 

0 

0 

1 

1 

0 

—t 

0 

—t 

1 

1 

0 

0 

—t 

—t 

—t 

—t 

0 

0 

1 

1 

—t 

0 

—t 

0 

1 

1 

0 

0 

0 

0 

1 

1 

—t 

0 

0 

—t 

—t 

—t 

0 

0 

1 

0 

—t 

0 

—t 

1 

0 

0 

0 

0 

1 

0 

—t 

—t 

0 

0 

0 

—t 

-t 

1 

1 

-t 

0 

—t 

0 

1 

1 

0 

0 

0 

0 

1 

1 

0 

—t 

—t 

0 

0 

0 

—t 

—t 

1 

J_ 

0 

—t 

0 

—t 

1 

J_ 

0 

0 

0 

0 

1 

_l_ 

—t 

0 

0 

—t 

-t' 

“ 0 

—t 

~0 

1 

-t' 

—t 

0 

~0 

1 

-t 

~ 0 

0 

-t 

1 

lf ' 

“ 0 

0 

~0 

0 

—t 

0 

—t 

1 

1 

—t 

—t 

0 

0 

1 

1 

0 

—t 

—t 

0 

1 

1 

0 

u 

0 

0 

-t 

0 

—t 

0 

1 

1 

0 

0 

—t 

-t 

1 

1 

0 

—t 

—t 

0 

1 

1 

0 

0 

u 

0 

0 

—t 

0 

—t 

1 

0 

0 

—t 

—t 

1 

—t 

0 

0 

—t 

1 

0 

0 

0 

u 


It is straightforward to identify the Ra'. 


0^x4 

-t{I + R4 

-t{l + Ry) 

—t{I + Rx) 

-t{I + R,) 

0^x4 

—t{I + Rx) 

-t{l + Ry) 

-t{l + Ry) 

—t{I + Rx) 

0^x4 

-t{I + Rz) 

—t{I + Rx) 

-t{l + Ry) 

-t{I + R4 

UI 


Given this form, it is natural to look to a basis. Working hrstly with [1,1,1,1]^, we look 

16-eigenvector based on one of the vectors of the for a solution to 


0^x4 

-t{I + Rz) 

-t{l + Ry) 

—t{I + Rx) 
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—t{I + Rx) 

-t{l + Ry) 

-t{I + Rz) 

UI 


-l 

T 



T 


1 



1 

a 

1 


a 

1 


_1_ 



_1_ 


T 



T 


1 



1 

b 

1 


b 

1 


1 



1 



= E X < 
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_1_ 



_1_ 
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1 



1 

d 

1 


d 

1 

- 

1 


\ 

1 
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where a, b, c and d are to be determined. Ac- 
counting for the action of each block on [1,1,1,!]^, 
which is Ra eigenvector with eigenvalue +1, we 
have the equivalent to the 4-dimensional problem 
Aa = {—E/2t)awith 


/0 11 1 \ 

/a\ 

10 1 1 


b 

110 1 

, a — 

c 

Vl 1 1 -U/2t) 

\dj 


Pollowing tedious calculation or the use of a com- 
puter, the characteristic polynomial of A is deter- 
mined to be 

x{x) = —U (l6t^ — 12t^x + x^) 

- ASt^ + m^x - 2At^x^ + x^. 


This confers four solutions. To obtain information 
of the ground state energy and its 1/U correction 
{U ^ t), scale X by a factor oil/U, as in 


where 


x-P3 + % 


P 2 >{x) = — (l6t^ — 12 t^x + x^) 
= —{x + 4:t){x — 

and 


P4{x) = + 64t^x - 24 ^ 2^2 + x^. 


When U = oo, Ps gives the least root x = —At. For 
hnite U ^ t, this is subject to a small correction ac- 
cording to 




■ X = —At + 


= -4t- 


U{x - 2t)2 
12^2 
TT’ 


X=—4:t 


which gives us a candidate for the ground state en- 
ergy. It may be verified that the other three vectors 
of the Ra basis in fact yield 4-by-4 matrices with the 
same characteristic polynomial. 


x{x) = {x — 2t)^(x + 2t)^ — Ux{x — 2t){x + 2t), 


which has least root x ^ —2t for large U. There- 
fore our first solution does describe the ground state, 
which is hence non-degenerate and of energy 


Eo = —At — 


12U 

TT 


We can go further to determine the corresponding 
eigenvector in the limit —)► oo. To do so, note that 


in order to satisfy the eigenvalue equation in this limit 
the bottom entry doia must be 0. This means we can 
‘knock out’ the final row and column of A and solve 
/0 1 1 \ 

1 0 15 = {-E/2t)a. 

Vl 1 0/ 

When E = — 4t, we obtain the eigenvector a = 
(1,1,1, 0)^ which corresponds to 

\Eo)=E\i} 

i=l 

in the original basis. This solution has important 
physical signihcance in relation to resonating valence 
bond (RVB) state proposed by R Anderson in 1987 
to account for high temperature superconductivity in 
the cuprates [1]. 


Conclusion 

Symmetries are an integral part of any mathemati- 
cal physicist s problem-solving toolkit. Here we have 
seen the symmetries of the tetrahedron partition a 16 
dimensional matrix H into 4-by-4 blocks, reducing 
the dimensionality of the problem by a factor of 4. 
For more complicated systems the reduction can be 
far greater and indeed may ofier the only way into an 
otherwise intractable problem. 
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Review of high temperature superconductmty and 
the resonating valence bond theory. 
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Smooth Piecewise 
Polynomials 

Josef Greilhuber 


When I was in my second year at university, my 
friend Sergey Yurkevich asked the following ques- 
tion: 

Is there a smooth hmction on R which is 
not a polynomial, but for any given point 
X G M, the sequence (x), n G N has 
only hnitely many nonzero terms? 

We spent quite some time on this problem, attempt- 
ing to construct examples as well as chasing ideas for 
proofs of the converse, e.g. via Baires theorem. We 
gained no ground there. But in the search for exam- 
ples, I found a peculiar hmction which almost (in the 
Lebesgue sense) answers the question - because it is 
a smooth piecewise polynomial. 


Proof. Notethatthe/cSatishes/c(l—x) = 1 —/c(x). 
We construct /c by a hxed-point iteration. 

For this, we dehne : C[0,1] —)■ (7[0,1] as 

Tf{x) = ^£ f{T{s)) ds, 

where T : [0,1] —)► [0,1] is given by 

{ 3x ^ ^ I 

1 l<x<l 

3 — 3x ^ < X. 


Construction 

This hmction turns up most naturally as a solution 
to a hmctional equation involving hrst derivatives, 
similar to the (somewhat) famous Pabius distribution 
hmction [ 1 ]. We hrst show that it exists and is unique, 
then show its piecewise polynomial properties 


Then if /c is a hxed point of T, it has the desired 
derivative on (0,1). 

Let U be the closed athne subspace of C[0,1] satisfy- 
ing /(0) = 0 and /(1 — x) = l — f{x) for all f ^U. 
Then U is invariant under T. 

We observe that f{s) ds = since 


Theorem 1. There exists a unique function fc G 
(7^(R) with /c(x) = 0 forx <= 0, fc{x) = 1 for 
x>l and 


f» 


'lf{dx) x<\ 

- 1 < r 

2 3 — 

^ |/(3 — 3x) I < X. 




s)) ds 


= 1 - 



ds. 


andhencethat f{T{s)) ds = |. 
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Now we observe that T commutes with reAection 
on [0,1]: 

= ^£/r(i-i) di 

-Ijjmo. 


\Tf{x)-:Fg{x)\ = - 


r {fT{s) - gT{s)) ds 
Jo 

r \fT{8) - gT{s)\ ds 
Jo 

/' \fT{s) - gT{s)\ ds 
Jo 

^Jo 

\j\f{s)-g{s)\ ds 


<3 
- 2 


<3 

2 70 


<\ sup \f{s)-g{s)\, 
^ S6[0,l] 


= 1-J^f{x), 
and hence TU C U. 


So is a contraction with A = By Banachs fixed 
point theorem, there is a unique hmction satisty- 
ing Tfc = fc- We obtain our desired hmction by 
extending this onto R in the obvious way. □ 


We now show is a contraction on U, and hence ob- 
tain a fixed point. We must show that for f,g^U 
there exists A such that 

sup \J^f{x) - Tg{x)\ < X sup \f{x)-g{x)\ 

xG[0,l] 1CG[0,1] 

But by symmetry, we only need to check for x < ^: 


Note that fc is differentiable at 0, since 

And //(0) = 0. By symmetry, //(1) = 0, and there- 
fore, / is continuously differentiable, and satishes all 
conditions imposed in the statement of the theorem. 
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A Pigure 1 : The function / and its first few derivatives (scaled). 


It now evident that / is in fact smooth, since arbi- 
trary derivatives can be computed piecewise from the 
hmctional equation, and /^(0) = /c(l) = so they 
agree wherever two pieces meet. 

The formula fc{x) = |/c(T(x)) allows us to com- 
pute by induction that 



n—1 


k=l 


wherever T^{x) is dilTerentiable around x for 
/c = l, 1. Note that {T^y{x) is piecewise 

constant. 


In other words, for n, /c G N with n > 1 and 
k < there exist a polynomial Pn.k degree 

at most n such that 


fc{x) =Pn,k{^) 


3/c + 1 3/c + 2 

gn ’ gn 


Notably, the set of points where / does not agree with 
a polynomial of hnite degree is exactly the cantor set, 
and thus on a dense open subset of R /c is piecewise 
polynomial. It can be given as a disjoint union of 
intervals In^k^ where n G N and 1 < k < 
together with (—oo, 0) and (1, oo). On every 
there exists a polynomial of degree exactly n, such 
that f\j- ^ = Pn,k^ justifying the titular remark that 
/c really is a smooth piecewise polynomial. 


References 


SinceT"(a:) = lfora: € (2|^, 2|^). the (n + l)«* 
derivative of /c vanishes on (^|^^, ^|^^)> implying 
that / agrees with its Taylor approximation there, 
which we denote by 
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The Gem of 
Combinatorks 

Trevor Cheung 


The associativity axiom of group theory is something 
so familiar to all of us that we seldom think about it. 

a ■ {b ■ c) = {a ■ b) ■ c 

It basically tells us that for three elements, the order 
in which we bracket them does not matter as long as 
they are not reordered. What about four elements, 
or in general n elements? It is not very dithcult to de- 
duce from the axiom using induction that all the ways 
of bracketing are equivalent. Here, we ask something 
of a “meta-question”: 

“How many ways are there to 
bracket (n + 1) elements?” 

Lets denote the answer to this as Cn, as it clearly de- 
pends on n. Although the above dehnition is most 
clear for n > 2, we interpret it also to mean 

Co = = 1. 

Next, we have 

(C2 = 2. 

To get a glimpse of how we can solve the question for 
general n, lets consider the case of n = 3. With a bit 
of thinking, we find C3 = 5 difierent ways of brack- 
eting: 

(ab)(cd) a(b(cd) a((bc)d) 

((ab)c)d (a(bc))d 

From here, we omit the binary operation symbol for 
brevity. 


Notice that we can put these bracketing styles into 
three categories, depending on the position of the 
outermost operation. With {ab){cd), the last oper- 
ation is between b and c; with a{b{cd) and a{{bc)d), 
it is between a and b; with {{ab)c)d and {a{bc))d, the 
outermost operation is between c and d. The reason 
why there are exactly two ways when the last oper- 
ation is between a and b or between c and d is that 
C 2 = 2: thats in how many ways we can bracket the 
remaining triple. 

Let s look at the case of n = 4 before we lay down the 
way to calculate Cn in general: 


Category 

Bracketing 

a ■ b 

a ■ {{bc){de)) 
a ■ {b{e{de))) 
a ■ {b{{cd)e)) 
a ■ {{{bc)d)e) 
a ■ {{b{cd))e) 

b ■ c 

{ab) ■ {c{de)) 
{ab) ■ {{cd)e) 

c ■ d 

{a{be)) ■ {de) 
{{ab)e) ■ {de) 

d ■ e 

{{ab){ed)) ■ e 
{a{b{ed))) ■ e 
{a{{be)d)) ■ e 
{{{ab)e)d) ■ e 
{{a{be))d) ■ e 


Again, the reason for having 5 ways of bracketing in 
the first and last category is Cs = 5. For the second 
and third category, where we break the five elements 
into a group of 2 and a group and 3, we have just one 
way of bracketing the 2 elements, but 2 ways for the 3 
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elements. So these categories each have Ci • C2 = 2 
elements. From this more involved example it is easy 
to see that, in general, 

Cn+l = CoCn + • • • + CnCo, 

or equivalently 


Cn+l — E CrCn-r 


( 1 ) 


r=0 


A closed expression 

To obtain a closed expression for Cn from (1), we can 
use the generating hmction of Lets dehne 


/(®) := y^,Cnx"-. 

n=0 


and consider [/(x)]^: 


(1 — 4x) 2 = 

00 

-E 


°° +H) 


2r-3 \ 
2 ) 


r=0 


= i + E(-i) 


2r 


{-Axy 

1 • 1 • 3 • 5 • • • (2r - 3) 


2 ^r\ 


( 4 *) 


(2r-2)! 
r!(r — 1)!"^ 

r=l ^ ^ 

1 - y\-^x = 


- 1 - 2 E 



(3) 


[/(*)]" = E 

n=0 

+ • • • + C„(7o) 

00 

= y^c7„+ix" 


n=0 


n=l 

00 

= -y2cnx^ 

X 

n=l 


- 1 


x[f{x)]'^ = f{x) - 1 


x[f{x)f - f{x) + 1 = 0 


( 2 ) 


We can apply the quadratic formula to (2) and obtain 


f{x) 


1 — \/\ — \x 
2 x 


as the other solution diverges when x H 0. To re- 
cover Cn, we will use the Taylor (or binomial) series 
of f{x): Note that 


And by simply comparing the coethcients, we have 




n + 1 


The Catalan numbers 

This number sequence is called the Catalan num- 
bers. Far from just a curious question related to the 
associatmty axiom of group theory, it is of vital im- 
portance to combinatorics because of its ubiquity in 
combinatorial problems. Here are two other, very 
much different-looking problems that also have Cn 
as the answer: 

(i) How many north-easterlypaths (paths thatgo north 
or east in each step)from (0,0) to (n, n) are there not 
crossing the line x = y? 

(ii) How many ways are there to triangulate a regular 
{n-\-2)-gon (i.e. slice thepolygon into n triangles using 
a selection of non-intersecting diagonals)? 

Question (i) can be illustrated with a diagram: In Fig- 
ure 1, the red line is x and the green path does 
not cross the red line line. Question (i) simply asks 
how many such green paths there are. 
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A Pigure 1: lllustration ofQuestion (i) 


It is possible to create a bijection between paths in this 
question and ways of bracketing (n + 1) elements. 
To obtain the corresponding path from a bracketing, 
we hrst move one step to the right, then traverse the 
bracketing (as a string of symbols), going to the right 
whenever we encounter an opening bracket and up- 
wards whenever we encounter a binary operation. It 
is possible to reverse this construction, so it is bijec- 
tive. 

For example, the path above corresponds to this way 
of bracketing 6 elements: (a • (6 • (c • d)) • (e • /). 


A Pigure 2: lllustration of Ouestion (ii) 

Pigure 2 shows several triangulations of a pentagon, 
as mentioned in Question (ii). It would again be pos- 
sible to construct a bijection, but we will instead ver- 
ify that the recurrence relation holds. Technically, 
there is no 2-gon, but Ci = 1 conforms to a trian- 
gle being already triangulated, while the recurrence 
relation can be explained as foIIows: For a general tri- 
angulation of an {n + 2)-gon, we choose one special 
edge and consider the triangle that it is adjacent to. 
This triangle splits the (n + 2)-gon into a (/c + l)-gon 
and a (n — /c + 2)-gon, each with independent and 
arbitrary triangulations. So the number of trian- 
gulations satishes 

n 

TrTn—r- 

r=0 


Tn+l - X] 


It has the same initial values and recurrence relation 
as the Catalan numbers, hence is identical for all n. 


Hankel transforms 

One of the most interesting properties of Catalan 
numbers is their Hankel transform. The Hankel 
transform is a transformation of an inhnite (integer) 
seguence ai, a^, as,... to another inhnite (integer) 
seguence /ii, /^2, /^3,, • • • by the foIIowing rule: 



/ai 

a2 

an \ 

det 

a2 

as 



\f^n 

an-\-i 

a2n-l/ 


For example, the Hankel transform of the Fibonacci 
numbers 1,1,2,3,5,... is the eventually constant 
seguence 1,1,0,0,0,.... WecaneasiIyprovethisfor 
n > 3 by observing that each row is a linear combi- 
nation of the previous two. 

Interestingly, the Hankel transform of Co, Ci, C2, •.. 
as well as Ci, (^2, Cs,... is the constant seguence 
1,1,1,.... There is a number of reasons why this is 
interesting: 

(a) It is not guaranteed that by shifting a se- 

quence, the Hankel transform stays the same (e.g. 
^^2, 0/4 ,... does not have the Hankel transform 

of a constant sequence (although it is actually a se- 
quence we all know and Iove!). 

(b) The fact that the determinants stay the same when 
we add more entries is pretty intriguing. 

(c) Although the Fibonacci sequence also hillils both 
(a) and (b) (except that it is only an almost constant 
sequence), this is clearly to be expected from our con- 
struction of the sequence; but for Catalan numbers, it 
is not obvious that its Hankel transform is a constant 
sequence even from the recurrence relation (which 
at hrst glance does seem to be the key to prove it, but 
actually is not). 

It is possible to use elementary combinatorics to ex- 
plain why the Hankel transform of the Catalan num- 
bers is a constant sequence, but this is too invoIved 
to be put in a short article. Needless to say, Cata- 
lan numbers are a real gem of combinatorics, soIving 
many seemingly unrelated combinatorial problems, 
as well as having interesting properties. 
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Doodles 

Jeremy Taylor 


A doodle, that is a curve that you can draw in the 
plane without lifting your pencil. 

The doodle divides the plane into faces. Notice that 
they can be coloured with two colours so that adja- 
cent faces have different colours. We will give three 
proofs that the faces of a doodle can always be 2- 
coloured. 


Euler Cycles 

Pirst we need to develop the right language to talk 
about doodles. A doodle can be formalized as a 
planar graph U, consisting of vertices connected by 
edges and embedded in the plane. And a doodle di- 
vides the plane into faces. 


But not all planar graphs correspond to a doodle. For 
example, the planar graph below is not a doodle. No- 
tice that the faces cannot be 2-coloured. 



Definition 1. An Euler cycle isapath that visits every 
edge exactly once. 

So a planar graph is a doodle if and only if it has an 
Euler cycle. Such graphs are characterized by the fol- 
lowing simple condition. 
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Proposition 1. A connected graph has an Euler cycle 
ifand only ifall vertices have even degree. 


Proof Suppose G has an Euler cycle 7. Then every 
time 7 enters a vertex it must also leave, so the de- 
gree is even. 

Now suppose that every vertex has even degree. We 
would like to construct an Euler cycle. Start at some 
vertex v and repeatedly follow unvisited vertices. 
Since all vertices have even degree we can only run 
out of vertices to follow after returning to v. So we 
have constructed a cycle. If it does not visit all edges, 
we can append additional loops until all edges are vis- 
ited. □ 


Two-Colouring 

Now we need to understand precisely the statement 
that the faces can be 2-coloured. 

Deiinition 2. Ehe dual of a planar graph is obtained 
by replacing faces with vertices and connecting those 
that are adjacent faces. We denote the dual of G by 
G\ 



So the faces of D can be 2-coloured if and only if the 
vertices of D* can be 2-coloured. 


Proposition 2. A graphs vertices can be 2-coloured if 
and only ifthere are no odd length cycles. 


Proof. We wish to colour the vertices black and white. 
Eirst we take some initial vertex v and colour it black. 


cO 



Then for any other vertex w choose a path 7 from v to 
w. Since the colours of vertices along the path must 
alternate, the colour of w is determined by length^y) 
mod 2. However we need to know that the result is 
independent of the path we choose. This happens if 
and only if there are no odd-length cycles. □ 


There are two themes to notice here. Eirst is par- 
ity. We are working mod 2. Second, there is an in- 
teraction between the local and the global. A two- 
colouring of a graph is determined locally. The con- 
dition that the graph has no odd-length cycles means 
that there are no global obstructions. 

Two-Colouring Doodles 

Theorem 1. The dual graph of any doodle D can be 
2-coloured. 


Proof. We have the following chain of reasoning. 

1. Disadoodle. 

2. D has an Euler cycle. 

3. All vertices of D have even degree. 

4. The dual D* has no odd length cycles. 

5. The dual D* can be 2-coloured. 

The only implication left to prove is 3 ^ 4. 

We say that a cycle in D* is simple if it goes around a 
single vertex. Suppose that all vertices of D have even 
degree. Then any simple cycle in D* has even length. 
















Below is a simple cycle around a vertex of degree 6 . 



Now suppose that we have a cycle 7 in D* that is not 
simple. 



Doodles on the Torus 

Now suppose we draw our doodles on the torus in 
stead of the plane. 



Then our cycle bounds a region containing a collec- 
tion of vertices, Vi.. .Vr. 



We break up the full cycle into simple cycles 71 ... 7r 
that bound the individual vertices. 

Then 

length 7 ^ = 2 / + length^, 

i 

where / is the number of interior dual edges (of some 
7i but not 7). 

Since length 7^ is even, we have that length 7 is even 
as desired. □ 


Consider the following counterexample. 



Gluing opposite sides of the square gives a doodle D 
that cannot be 2 -coloured. 

So we need to see where the argument goes wrong. 




In the counterexample there is an odd length cycle 
in the D^. The problem is that our odd length cy- 
cle goes around a hole torus so it does not bound a 
region. Thus the argument in Theorem 1 fails. 





























Winding Number 

DeAnition 3. For any point p G that does not lie 
on the doodle D, the winding numher I{yD^p) is the 
numher ojtimes that the curve 7^ wraps around p. 



Notice that I^yn^p) is constant on the faces. Fur- 
thermore the winding number of adjacent faces dif- 
fers by ± 1. This gives a 2-colouring of the faces where 
the colour is determined by the parity of the winding 
number. 

Proposition 3. The winding numher differs hy±lon 
adjacent faces. 

Proof. 



We wish to prove that the winding numbers, 
I{jD,Pi) and/( 715 ,^ 2 ) diiFerby ±1. 



To do this we consider adding a path 6 around Pi that 
cancels out’ part of the edge between Pi and P^. We 
have 


I{7D + Pl) = I{lD. Pl) + m Pl) 

= 7(7d,A)±1. 



But after adding (5 to 7 ^ we have that Pi and P 2 lie 
in the same connected component. So 

I{1dPS,Pi) = I{idPS,P2) 

= I{lD,P2)+I{S^P2) = I{lD,P2)- 

Hence I{id: + 2 ) = Hid, + 1 ) ± 1 as desired. □ 


Cycles and Boundaries 

Consider a doodle D. 

Let V be the set of vertices, E the set of edges, and F 
the set of faces. 

We dehne 

C7o = 01"iZ/2, Ci = 077*2/2, 

and C72 = 0FiZ/2, 

the free Z /2-modules generated by the set of vertices, 
edges, and faces respectively. Note that each element 
of Co, Ci, C 2 corresponds to a subset of V, E, F re- 
spectively. 

Then there are ‘boundary’ maps (Z/2-module homo- 
morphisms) 

772 Ci > Co. 

The map sends a set of faces to the set of edges on 

its boundary. 



















The map di sends a set of edges to the set of vertices 
on its boundary. 


An Example 


29 



C, 


We have 


Consider the following planar doodle where we have 
labeled the vertices, edges, and faces. 


E, 



di o d 2 = 0 

or equivalently 

im d^ C ker di , 

because the boundary of a boundary is zero. 


For example, we have di{Ei) = Vi + V 2 are the two 
vertices on the boundary of Ei. 

Similarly ^ 2 (^ 2 ) = E^ + Ei^ + E^, + E 20 are the 
four edges on the boundary on E^. 

As expected these form a cycle, so di o d^^E^) = 
+ ^15 + ^3 + ^20) = 0 - 


The elements of ker di are called cycles because th( 
correspond to (one or multiple) closed loops on tl 
graph. 



The elements of im d^ are called boundaries becau 
they correspond to loops that bound some collectic 
of faces. 




Another example of a cycle is B = Ei -\- E^ -\- E^ + 

Ell,Ei2 + T^13 -\- Eq -\- Ei8 G Ci. 


cyc\e) 


We have di{B) = 0 so B is indeed a cycle. 

But B is also a boundary since B = d^^Ei + +9 + 
+8 + +"11)- 



















30 


Two-Colouring Doodles 

Suppose that D issi doodle in the plane. Then since 
the plane has no holes, every cycle bounds some set 
of faces. So we have equality im = ker di. 

Theorem 2. Jhe dual graph ofa planar doodle can be 
2-coloured. 


D but only on the topological space X in which the 
graph is embedded. Given a space X, we can thus de- 
fine iTi(X, Z/2) which encodes information about 
the 1 -dimensional holes in X. 

We said before that since the plane has no holes 
iJi(R^,Z/ 2 ) = { 0 } istrivial. 

On the torus, 


Proof Consider the sum of all edges 

B = Y,EieC2. 

We have di{B) = 0 since every vertex has even de- 
gree. Therefore since ker di = im d^ it follows that 
B G imd^. 

Let 

B = d 2 {A) for some subset A = OiPi G C^. 
Then gives a 2 -colouring of the faces. □ 


iLi(r2,Z/2) ^ Z/2 0Z/2 

which corresponds to the two holes of a torus. 



Spaces with Holes 

We have said that im d^ = ker di for any graph D 
embedded on a topological surface without holes. On 
a general space we do not always have equality. A cy- 
cle will fail to be a boundary when it goes around a 
‘hole. 

To be more precise, recall that the cycles, ker di, cor- 
respond to loops. But the boundaries, im are the 
cycles that don t actually go around a hole so we quo- 
tient out by them. 

Thus we define the first homology module 
Hi{D, Z/2) = ker di/ imd^. 

As long as the faces of D are reasonably behaved, 
iLi(ii,Z/2) doesnt actually depend on the graph 


For instance, 7 i = {Ei, Es} G ker di \im is a 
cycle. But it is not the boundary since it goes around 
a hole. 


J 





1 

1 

' 

J 




u 



t 



Also 72 = {E^, E^,Eq} G ker di \ im d^ is a cycle 
but not a boundary. 

Thus [ 71 ] and [ 72 ] represent nonzero elements in 

Hi(r2,z/2). 

Purthermore they are generators, 

Hi{T\Zl2) = ([7i],[72]) ^Z/2©Z/2. 



Other Coefficients 

We can also do homology with coethcients in other 
rings. More generally for a ring R, we can define 

CQ = ^ViR, Ci = ^EiR, 

and C 2 = 0Fii?. 
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This just means that every element of Cq is of the 
form a sum of vertices with coehicients in 

R. In the case R = elements of Ci correspond 
to subsets, so we recover our old dehnition. We only 
work with coethcients mod 2 so we do not have to 
worry about orienting faces. 

In the general case, orienting the faces and vertices is 
necessary to get boundary maps 

C 2 > Ci > Co with imd2 C kerdi. 

Again we can dehne 

Hi {G, R) = ker di / imd^. 


Two More Colouring Results 

Theorem 3. Let Gbe a cubic planar graph. 

Jhen thejaces can be L-coloured ijand only ijthe edges 
can be 3-coloured. 

Prooj We will consider colourings of the faces, ver- 
tices and edges by the four elements ofZ/20Z/2. 

Then colourings of the faces correspond to elements 
of C 2 . And colourings of the edges correspond to el- 
ements of Ci. (Here we do not yet require that adja- 
cent faces or edges have different colours.) 

Let A G C2 be a colouring of the faces and B ^ Ci 
be a colouring of the edges. 

1. Adjacent faces have different colours in A if 
and only if d^^A) is a colouring of the edges 
by nonzero elements ofZ/20Z/2. 

2 . Adjacent edges of B have different nonzero 
colours if and only if di {B) = 0. 

Suppose that A is a 4-coIouring of the faces. Then by 
I we have that d^ {A) colours the edges with the three 
nonzero colours. But di o d^^A) = 0. So by 2 we 
have that d^^A) assigns distinct colours to adjacent 
edges. 

For the other direction we use the fact that the plane 
has no holes so Hi{G^ Ij/I 0 Z/2) = {0}. 

Suppose that B is a 3-coIouring of the edges. Then 2 
implies that B G ker di is a cycle. And Hi (G, Z/2 0 
Z/2) = {0} means that B = d^^A) is the bound- 
ary of some A G C2. By I, adjacent faces of A have 
different colours. □ 


Theorem 4. Let Gbea triangulation oja simply con- 
nected surjace. Then the vertices can be 3-coloured ij 
and only ijC is a doodle. 

Prooj. 

(^) First suppose that G is not a doodle. Then G has 
a vertex Vi of odd degree. 



But Vi and its adjacent vertices cannot be 3-coIoured. 

(<^) Now suppose that (T is a doodle so the faces can 
be 2-coIoured. Orient edges based on the colouring 
of the faces (so that the the black face is always on the 
right). 



Now we do cohomology with coellicients in Z/3. A 
l-cocycle is a function on edges. Such a l-cocycle in- 
duces a functions on faces by summing the values on 
the oriented boundary edges. For example, consider 
the sum ofthe coedges /3 :=Y/ which is the func- 
tion that takes the value I on all edges. Its cobound- 
ary vanishes modulo three since the faces are oriented 
triangles. 

Since the space is simply connected, our l-cocycle 
[3 must be the coboundary of a 0-cocycIe, Y — 
^(a) for some ol = Y ^ means there 

is a function a on the vertices such that the difference 
between its values on adjacent vertices is I. There- 
fore the values G Z/3 of the 0-cocycIe give a three 
colouring of the vertices. □ 

Using the dual of the graph, the foIIowing result can 
also be obtained. 

Theorem 5. Let G be a planar cubic graph. Then the 
jaces can be 3-coloured ijand only ijthe vertices ojG 
can be 2-coloured. 







THE QUANTITATIVE TRADER PR06RAMMES 


SIG's Quantitative Trader Internship and Graduate Progrannnnes can 
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DURING THE INTERNSHIP^ YOU WILU 

Work directly with our Quantitative Traders, learning 
how they make trading decisions and assess uncertainty 
in the markets 

Learn how to make rational decisions by participating 
in mock trading and playing strategic games. This is a 
great way to explore how our traders go through similar 
thought processes while evaluating the expected value 
of a given trade and determining how to price risk. 

Find out how we use game theory to hone our skills and 
apply them to complex capital markets 

Participate in daily education where you will learn option 
theory and how to make positive expectancy decisions in 
the financial markets 

Work on a specific trading desk and complete a project to 
help your desk tackle unsolved problems (like analysing 
trade data, improving efficiency, or creating new tools 
for traders) 

Develop new skills using the technology tools our 
Quantitative Traders use, including Python, SQL, and VBA 


DURING THE GRAD PRDGRAMME, YDU WILL 

Work as part of a trading desk 

Complement your on-the-desk learning with 
classroom education 

Expand your knowledge of decision science, game 
theory, and capital markets 

Participate in mock trading to reinforce what you are 
learning on the trading desk 
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A Four-Player Game of 
lntransitive Dice 

Jared T Collins and Isabel E. Harris 


We produce a set of nineteen intransitive dice satisfy- 
ing that for any subset of three dice there exists a die 
which will roll higher than each of them with proba- 
bility greater than one half. Moreover we do so using 
nine-sided dice which is currently thought to be the 
least possible and allows this game to be physically 
created. 


PoMtely Cheating 

Transitivity is one of our most prevalent mathemati- 
cal properties: A > B and B > C then A > C. 
Unfortunately transitivity is a terrible property for 
games. Imagine a transitive variant of Paper-Rock- 
Scissors. Paper beats Rock and Rockbeats Scissors so 
Paper beats Scissors too. Paper is king in this game, 
and you would be foolish to pick either of the other 
items. 

No, we want transitivity far away from the mechanics 
of our games. Not to say we want it completely re- 
moved: scores should still rank players when needed. 
But you will be hard-pressed to find people willing to 
play a game in which one specihc move guarantees a 
win. So instead of guaranteeing a win, can we just tip 
the odds in our favour a little bit? 

Consider the following game. We both roll a standard 
six-sided die and the higher number wins. Seems 
pretty boring, but I bet you have played it before to 
decide who goes first in some other game. How about 
we play a more complicated game? Instead of stan- 
dard six-sided dice, we will use a special set of dice. 
The numbers on their sides are the following: 



A casual look at these dice suggests another boring 
game, but at least there will not be any ties. The dice 
all have an expected value of 4 and the dice often re- 
peat. But if it were a boring game we would not be 
talking about it. A closer look shows that die 3 beats 
die 2 five out of nine times, die 2 beats die 1 five out of 
nine times and die 1 beats die 3 five out of nine times. 
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The relative strengths can be expressed as: 



So we have an intransitive system much like Paper- 
Rock-Scissors but with some random chance thrown 
in to make us look less like cheaters. But, to make our 
scheme complete we will need to ensure we always 
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pick our die last. That is why we are polite cheaters: 
always allowing our opponent to pick their die tirst. 

Doublethe Fun 

There are several ways we can improve the game we 
have just described. Our winning probability of 0.5, 
while better than 0.5, could be higher and there are 
dice sets that achieve higher probabilities, although 
it is impossible to surpass 0.75 [6]. We could try to 
use fewer dice or reduce the number of sides on each 
die, but three dice is the best we can do in a two- 
player game and our dice are already effectively three- 
sided. Keep these three kinds of improvements in 
mind: higher probability, fewer dice and fewer sides. 

We, however, are concerned with a different kind of 
improvement. Rather than play with just one friend, 
could we play with more and cheat two people at 
once? It is not quite the same as playing a single 
person twice, since the probabilities against two op- 
ponents are not independent. But it should be clear 
that our set of three dice is not up to the task. Can 
we find a set of dice such that for any pair of dice 
there is another die in the set that beats each of them 
with probability greater than one half? This is a spe- 


cial case of the Schiitte Tournament Problem [4], 
which asks us to find a complete directed graph with 
n vertices such that for any subset of k vertices there 
is at least one vertex with edges to each of the k ver- 
tices and to do so with the smallest value of n. In our 
discussion the direction of the edge is equivalent to 
one die beating another with probability greater than 
one half. 

Such a set of dice exists. In 1992 Oskar van Deven- 
ter [7] described a set of seven three-sided dice that 
realize a Schiitte Tournament with k = 2. The dice 
beat one another with a probability of 0.5, just like 
our first example. 

Going Bigger 

Winning against two friends at a time is great. But we 
can do better. Inspired by Van Deventer s dice and 
with the knowledge of how many dice are needed for 
the Schiitte Tournament with k = 3 we set to work 
trying to create a set of dice with a small enough num- 
ber of sides so that we could physically create them. 
At the time we believed the best set had dice with 342 
sides [5], not something you would be able to hold in 
your hands. We got to work and made a set with 9 
sides. 



◄ Pigure 1 : Van Deventer's dice set 
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A Figure2: Oursetof 17dice 


▼ Pigure 3: Die 10 vs. Die 14 Pigure 4: A general match T 
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Of course, physically making a die with 9 sides is a bit 
of a challenge. So much like how Van Deventer dis- 
guised a set of three-sided dice by duplicating each 
side to create a six-sided die, we will need a trick to 
make our dice more standard. The easiest way is for 
us to embed each nine-sided die into a ten-sided die 
by letting the tenth side be blank or giving it a sym- 
bol that instructs the player to re-roll the die. So here 
they are: Pigure 2 shows the set we found. 

As you can imagine, in a complete digraph with nine- 
teen vertices it gets a little dithcult to see the direction 
of each arrow. So we have simplihed the picture to 
focus on one die and the dice that it beats. If you are 
curious to know what dice any specihc die beats you 
can determine that utilizing the rotational symmetry 
of the graph or by utilizing the quadratic residues of 
19. The probability associated with one of these dice 
beating another is l^. That is not quite as strong as 
the two- and three-player games, but it is still in your 
favour. Pigure 3 shows the comparisons between Die 
10 and Die 14. Pigure 5 shows a more general picture 
for any pair of these nine-sided dice. In each com- 
parison the main diagonal will split into 5 wins and 4 
losses, though the pattern of that split will vary. 

Can We Do Better? 

Think back to our three means of improving the 
game: higher probability, fewer dice and fewer sides. 
Add to those our fourth way to improve the game: 
more players. Higher probability is an area where 
there is plenty of room for improvement. || is very 
close to ^ and we know that the upper bound on pos- 
sible probabilities is |. Note that if we were willing to 
relax our condition of uniform probabilities this gets 
far easier to achieve. 

Pewer dice is an area that we cannot improve without 
adjusting the rules of the game. It is known that the 
Schiitte Tournament problem with k = 3 requires 
n > 19. However James Grimes [7] has had success 
in using fewer distinct dice in the three player game. 
His method does require you to have two of each die 
and the rules of the game are a little different but a 


similar technique could lower the number of unique 
dice needed for the four player game. 

Pewer sides is for us the most compelling avenue 
of improvement. Independently of our work Dezs6 
Bednay and Sandor Boz6ki [2], [3] also have created a 
set of nine-sided dice. They are working on the much 
larger problem of realizing any digraph with a set of 
intransitive dice. Levi Angel and Matt Davis [1] also 
have accomplished some great work in creating in- 
transitive dice sets for digraphs. This area is most in- 
teresting to us because any decrease in sides would 
necessarily bring an increase in probability. 

Pinally the issue of more players is interesting but is 
tied up in the Schiitte Tournament Problem [4]. All 
of the games described in this paper are examples of 
Paley tournaments [3]. If we believe that the optimal 
Schiitte tournaments are all Paley tournaments then a 
five player game would require 67 dice, but there may 
be a way to do it with fewer dice if we do not require 
our set to form a Paley tournament. 
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Braids, Knots and 
Grids 

Jessica R Wang 


To a large degree, the beauty of mathematics lies in 
the deep connections between seemingly ditferent 
helds. For millennia, knots have been used as a tool 
for fastening and nautical devices. Today knot theory 
is known as a signihcant held in topology. However, 
it was not until the 1860s when Lord Kelvin, a math- 
ematical physicist, hrst discovered a theoretical use 
of knots. He modelled atoms as knots in the aether, 
which led to the birth of topology after Peter Taits 
work of classifying knots [ 1 ]. Later, Emil Artin gave 
an explicit introduction to braid groups in 1925. I 
have drawn three diagrams for the cover picture, but 
they all can be seen as representations of the same 
thing, the left-hand trefoil knot. 


Braids 

In simple terms, a braid is a collection of strands at- 
tached to two opposite ends such that each strand 
never turns back towards its starting point, and the 
strings never penetrate one another. Eormally, in R^, 
let Ai = (i, 0 ,1), Bi = (i, 0 ,0) for 2 G {1, • • •, n). 
An n-braid is a conhguration of n strands, di, each 
running from Ai to some Bj, such that: (i). The 
strands do not intersect one another; (ii). Every level 
plane Pg = {(x, y, s)} parallel to the x^-plane inter- 
sects each strand at one and only one point. Eor ex- 
ample, in the second hgure below, the green strand 
turns upwards at one location, therefore it is not a 
braid. 



A Pigure 1 

A 3-braid 


A Figure2 

A non-braid 


The mathematical dehnition of a knot is a closed 
string in which does not self-intersect, and a link 
is a collection of knots. Two knots are equivalent (~) 
if one can be transformed to the other by a hnite se- 
quence of Reidemeister moves. The four types of Rei- 
demeister moves (Type 0,1, II and III) are shown on 
the right. Similarly, equivalent braids are ones which 
can be deformed into one another without cutting 
any strands. 



A Pigure 3: The three Reidemeister moves 
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A Pigure 4: The L-nnove on a figure-8 knot 


For simplicity here, lets refer both knots and links 
as knots. A braid [3 can be easily transformed into a 
knot by connecting its corresponding ends, forming 
a hraid closure / closed hraid 


can produce ditTerent braids. For example, applying 
the Type I Reidemeister move to downward -tlowing 
strands allows us to add more strands to the resulting 
braid: 



The reverse process, turning a knot into a braid, is 
less straight-Torward. But James Waddell Alexander 
II told us that it is always possible [2]: 


Theorem 1. Alexanders Theorem Every knot or link 
can he transjormed into a closed hraid. 


Intuitively, to prove Alexander s Tbeorem, one would 
move all the crossings of a knot to one side in an at- 
tempt to construct a braid closure. However, some 
strands might be oriented upwards. Tbe L-move cor- 
rects this problem by cutting each upward-running 
strand into two downward-running strands: 



Figure 4 shows this move applied to a figure-8 knot. 

Since a knot can be presented in different ways up 
to equivalence, applying L-moves to the same knot 



You could connect the corresponding ends of the 
braid in (5) to see that it is equivalent to the knot 
in (1). Note that each continuous section of the 
knot (between under-crossings) which points up- 
wards needs an L-move, and the number of L-moves 
required gives the number of strands in the result- 
ing braid. But how do we know that the L-moves are 
valid, i.e. the closure of the resulting braid is always 
equivalent to the original knot? We can use grid dia- 
grams to prove the vahdity. 


Grid Diagrams 

Knots are generally presented in a curvy structure. 
This nature means they are hard for a computer to di- 
gest. One way to increase computational ethciency is 
to represent a knot by a grid diagram, which only has 
vertical and horizontal strands, and two grid moves 
instead of the Reidemeister moves in knot theory. 
Grid diagrams are hmdamental in dehning grid ho- 
mology [3]. Agrid diagram is an n x n grid with ex- 
actly one o and one x in each row and column. It can 
be transTormed into a knot by connecting the o and 
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X in each row or column, with vertical lines always 
over horizontal ones. As for the reverse process, we 
hrst draw the knot so that the strands are horizontal 
and vertical, then arrange the vertical sections to be 
on top. For example, we have a grid diagram of the 
left-hand trefoil knot (Can you work out how many 
5x5 grid diagrams there are for the left-hand trefoil? 
Any generalisations?): 
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The grid diagrams are related by two grid moves, 
namely (de)-stabilisation and communication: Stabil- 
isation is the addition of a kink in any direction, while 
destabilisation is the removal. This preserves the as- 
sociated knot but changes the grid size. Below shows 
an example of adding a kink to the top-right of x. 


|x-| - ^o| 1 



X- 

=F^ 


-o 



Theorem, we introduce the mid-grid move [4], which 
resembles the L-move in knot diagrams. It breaks 
each upward strand into two downward strands, and 
creates an upward strand on the right edge of the grid. 
Each move increases the grid size by 2. In this way, a 
braid closure can be formed. 
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Proof of Alexander's Theorem 

The mid-grid move preserves the knot type as it can 
be expressed as a sequence of (de)-stabilisation and 
admissible communication moves. Consider an up- 
ward strand at the ith column in an n x n grid, with 
X and o on the jth and kth row respectively: 


Communication is the interchange of two consecu- 
tive rows or two columns. This preserves the grid size 
but can change the associated knot. A communica- 
tion is called admissible if the associated knot is also 
preserved. 
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If we dehne Q = grid diagrams, JC = 
equivalence classes of knots, then the map Q ^ JC 
is not injective, but we have a bijection 

JC i — Q! {(de)-stabilisation, adm. comm.). 

Now, back to our objective of proving Alexanders 


Step i: Add kink to the top-right of x; 

Step ii: Communicate (n — j) times to put o — x at 
the top; 

Step iii: Similar to i and ii, add kink to the bottom- 
right of o then communicate {k — 1) times to put 
X — o at the bottom; 
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Step iv: Communicate (n — i) times to put the (i + 
2)th column to the rightmost of the grid. Given a 
knot K, by the bijection we can transform it to a grid 
diagram, which then can be transformed into a closed 
braid by the mid-grid move. □ 


The Alexander Polynomial 

We have now proved Alexander s Theorem. Unsur- 
prisingly, our topologist Alexander also discovered 
the famous Alexander polynomial for knots. It is a 
knot invariant, meaning that if two knots have dilfer- 
ent Alexander polynomials, then they are not equiv- 
alent. But two dilTerent knots may have the same 
Alexander polynomial. 

An easy way of computing the Alexander polynomial 
^L(t) of a knot L is to hrst compute the Alexander- 
Conway polynomial V^(t)^ which is a modihed ver- 
sion of the Alexander polynomial found by John 
Conway [5]. It is dehned by Conways skein relation: 

^ unknot(t) 1 ? 

where the unknot denotes the simple closed loop ‘O’; 


Vl+ — Vl_ = t • Vlo, 

where L is any knot diagram and L^, L_, Lo are 
dehned by L^ = L; L_ for the knot diagram ob- 
tained by changing the crossing to a negative; Lo for 
the knot diagram obtained by eliminating the cross- 
ing such that the orientation is preserved. 

/\ /\ 

i_|_ ig i_ 


The Alexander polynomial Al (t) can then simply be 
calculated by the relation 

AL(f')= VL(f-/'). 

Try it with the left-hand trefoil knot, you should get 

Aleft-hand trefoil(f) t 1 + t 

Note that although the Alexander-Conway polyno- 
mial is considered as a simplihcation of the original 
Alexander polynomial, the latter has given signihcant 
results in knot Ploer homology [3]. 


The Alexander Polynomial using 
Grids 

As said before, grid diagrams make calculations more 
computationally possible, and this includes comput- 
ing the Alexander polynomial. The hrst step is to 
construct the Minesweeper matrix, which connects 
Seijert diagrams with grid diagrams. As we work 
through the process, it will become clear that its 
‘Minesweeper’ name comes from its unique struc- 
ture, which looks rather like the well-known com- 
puter game. 

A Seifert diagram of a knot is created by eliminating 
its crossings, which is done via connecting the un- 
derpass with the overpass such that the orientation is 
preserved: 



( 


This process creates closed loops; we’ll call them 
Seifert circles. For our left-hand trefoil, we have: 



Now, for each Seifert diagram with n Seifert circles, 
namely Ci,..., Cn, dehne the winduj(Ci) of Ci (i G 
{!,..., n}) as follows: If Ci is not in the interior of 
another Seifert circle Cj (i j), then if Ci has a 
clockwise orientation, uj(Ci) = 1; if \ has an an- 
ticlockwise orientation, cj((7i) = — 1. 

If however \ is in the interior of another Seifert cir- 
cle Cj (i j), then if Ci has a clockwise orienta- 
tion, uj(Ci) = uj(Cj) -\- 1; Ci has an anticlock- 
wise orientation, uj(Ci) = uj(Cj) — 1 . So if we take 
our previous diagram and label the bigger circle Ci, 
the smaller one C 2 , then we have uj(Ci) = 1 and 
^(C^) = 1 + uj(Ci) = 2 . 

Call the (i,j)th lattice point the point where the ith 
row intersects with the jth column, counting from 
the top-left lattice point. For each lattice point in the 
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Pigure 5: Knots, Grids, and Braids ▼ 
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grid, we can label it with the wind of the Seifert circle 
it is located inside. If it is not in any Seifert circles, 
then we simply write 0. Each of these are called the 
wind ofa lattice point, written as cj(2, j). Eollowing 
our left-hand trefoil example, we have 



Einally, we have set up the foundation of construct- 
ing the Minesweeper matrix! Eormally, if we let an 
n X n grid diagram G represent an oriented knot di- 
agram, then the Minesweeper Matrix M{G) of G is an 
n X n matrix {gij) constructed as follows: Eor each 
i ^ n — 1 and j ^ n — 1, we have gij = ^t all 

other points, g^j = 1. 

Basically, M(G) is a square matrix with t taken to the 
power of the individual winds as its elements, then 
being added an additional column and row of Is. 


So denoting the grid diagram of our trefoil knot as 
M{G) we get the following: 


M{G) 


/1 t t t 1\ 

t t^^ t^^ t 1 

t t^^ t 1 1 

t t 1 1 1 

Vi 1 1 1 1/ 


Knowing how to create the Minesweeper matrix, we 
are one step away from calculating the Alexander 
polynomial. In order to specify the locations of the o s 
and the x s, we use a systematic approach by letting 
O be the set of unit squares containing a o, and X be 
the set of unit squares containing a x. Then, counting 
from bottom to top and from left to right, suppose a o 
is in the intersection between the ith column and the 
jith row. Let ao be the permutation which maps i to 
ji, such that (70 = (ji j 2 • • • jn)- crx is dehned in 
a similar fashion. Moreover, let e(0) denote the sign 
of permutation which connects cro with the permu- 
tation {nn — 1 ... 1). Lastly, consider our previous 
grid diagram which has the associated winds. Notice 
that each o or x is contained in a unit square with 
wind(s) on its four corners. Let a(G) be the sum of 
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all the winds on the four corners of each o and x, di- 
vided by — 8. However, remember that our mission 
is to make things more computationally effective, one 
way is to let, say, Wi be the set of winds which are in 
the corner of exactly one o or x, be the set of 
winds which are in the corner of exactly two o s or 
X s, and so on. Using this method, we obtain 

a(G) =(a:- ^(*’-^') 

k=l \ uj{i,j)eWk 


of exactly one o or x in red; and the wind which are 
in the corner of exactly two o or x in blue: 
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Using our dehnitions above, we get our hnal result. 
The Alexander polynomial of a knot K is given by 

= e(G) • (r 5 - • det (M(G)) • 

where G is a corresponding grid diagram of K [3]. 

One way of proving it is to show /S.K(t) is in- 
variant under admissible communications and (de)- 
stabilisations. One would also need to prove that it 
satishes Conway s skein relation. 


From this we obtain 

a(G) = — \ ■ numbers in red 

— - numbers in blue 

A 4-^ 


= -F(6 + 2-9) 
= -3. 


The left-hand trefoil 

To summarise, we have proved Alexander s Theorem, 
which connects braids with knot theory, by using 
grid diagrams. Then we showed an alternative way 
of computing the Alexander polynomial, of which it 
would be natural to show an example: As you would 
have guessed, we are going to compute the Alexan- 
der polynomial of the left-hand trefoil. We need to 
determine n, e(G), det(M(G)), and a(G): 

n = 5, since we are working on a 5 x 5 grid diagram. 

e(G): From the leftmost column to the rightmost col- 
umn, the o s occupy the Ist, 2nd, 3rd, 4th and 5th 
rows respectively. So we have cro = (1 2 3 4 5). 
To obtain the sign of permutation which connects 
(70 with (n n — 1 ... 1) = (5 4 3 2 1), notice 
that the two are equal if we swap 1 with 5 and 4 with 
2. In other words, we would need to find the sign 
of (1 5) (4 2). Since this is a permutation with two 
transpositions, we have e(G) = 1. 

det(M(G)): We have already computed M(G) ear- 
lier, and 

det(M(G)) = — 5t^ + llt^ 

- Ut^ + llt^ - 5t + 1. 

a(G): We combine the grid diagram with the winds 
of its lattice points to give the below figure. For con- 
venience, we label the winds which are in the corner 


Thus, we have our desired polynomial: 

^left-hand trefoil(f) ~ 

= e(G) • (rl - i5)i-” • det (M(G)) • 

, - 14^3 + 11^2 -5t+l _o 

= 1-j-i- • t 3 

(t-2 -t2)4 

= (k -t^ + 1 ^) • t-3 
= t- l+1-^, 

which matches the result from our earlier method us- 
ing Conway s skein relation. 
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Topological Data 
Analysis 

Simon Heuveline 


Mathematics is the art ofgiving the same name to dif- 
ferent things. 

— Henri Poincarĕ 


Algebraic Topology 

This section gives a very informal introduction to the 
subject of algebraic topology. The goal of algebraic 
topology is to classify spaces by their global struc- 
ture. For this aim, we have to dehne what it means 
for two spaces to have the same global structure. We 
will focus on the very intuitive notion of homeomor- 
phic spaces. We say that two spaces Mi and are 
homeomorphic (written Mi = M^) if they can be 
continuously deformed into each other, i.e. there is a 
continuous bijection <f : Mi ^ M^ that has a con- 
tinuous inverse. The well-known standard example is 
that a doughnut can be continuously deformed into a 
coffee mug and hence they are homeomorphic. This 
deformation can be seen in Pigure 1. 

In a certain sense, it is an easy task to prove that 
two spaces are homeomorphic since we ’just’ have to 
write down a homeomorphism between them. It is a 
much harder task to prove that two spaces are non- 
homeomorphic since we would have to ’try all possi- 
ble homeomorphisms’ which is, of course, a nonsen- 
sical idea. Algebraic topology gives a way to tackle 
the latter task in a conceptually new way: We assign 
an algebraic object (e.g. a number, group, ring, etc.) 
to each space. Moreover, we demand that homeo- 
morphic spaces get assigned with the same (at least 
up to some notion of isomorphism) algebraic object. 


Such an assignment is called a (homeomorphism- 
)invariant. If one has properly dehned such an in- 
variant, one can compute the invariants of two given 
spaces and compare them. If they are not the same, 
we have proved that the two spaces are not homeo- 
morphic. Observations like this are the starting point 
of classihcation results for certain spaces. To give an 
intuitive example of such an invariant, the ’the genus 
of a surface’ we can consider the assignment 

jdi : {2-dimensional surfaces} —)■ N 

M ^ #{Holes of M}. 

If we can continuously deform a surface into an- 
other one (think about squeezing and stretching it, 
as in Pigure 1), we will not change the number of 


▼ Pigure 1: A homeonnorphisnn fronn a 
doughnut to a coffee mug. 
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holes, which makes a homeomorphism invariant. 
In this way, we could prove that a doughnut is not 
homeomorphic to a pretzel (see Pigure 2 ), because 
/3i(Pretzel) = 87 ^! = /3i (Doughnut). 

This is, of course, not a very interesting result. How- 
ever, once we formalise the steps of the heuristic dis- 
cussion above, we can also dehne a notion of holes 
and higher dimensional holes of higher dimensional 
spaces for which we have a priori no geometric intu- 
ition, so that the question of whether two spaces are 
homeomorphic may be highly nontrivial. This (and 
much more) is done in the framework of so-called 
homology and cohomology theories. The important 
invariants in homology, the Tomology groups’ (for 
us these will be in fact Q-vector spaces) Hk{X) can 
be formally dehned for a general topological space 
X. However, for us, it will be enough to know 
that l3k{X) = dimQ{Hk{X)), which is called the 
k-th Betti-number of X, counts the number of ’/c- 
dimensional holes of X’ in a certain sense and gener- 
alises the above intuitive notion of holes. An exam- 
ple of a ’one-dimensional hole’ would be the hole in a 
circle, while an example of a ’two-dimensional hole’ 
would be the hole in a 2 -dimensional sphere (think 
about a Kinder egg, with the ’hole’ being the interior 
of the egg). 

Even though algebraic topology is a very deep 
mathematical theory and applies to several inner- 
mathematical disciplines, it seems rather abstract and 
not suitable to be applied to ’real-world problems’ 
However, with time topology started to become a 


topic of interest in theoretical physics. There, it ap- 
pears on a hmdamental level in Quantum held the- 
ory and String theory but also in statistical physics. 
Eor instance, the genus (as discussed above) of a 
string worldsheet gives us the order in a perturba- 
tive expansions of string interaction amplitudes in 
analogy to the number of loops in the context of or- 
dinary Eeynman diagrams. Topology is actually a 
key discipline in the interplay between mathematics 
and physics since physical intuition turned out to be 
usehil for proving mathematical results about certain 
topological spaces (e.g. many results about 3 and 4 di- 
mensional manifolds or the local Atiyah-Singer In- 
dex Theorem). A rather new and quite surprising de- 
velopment on which I want to focus in this article is 
the application of algebraic topology to computer sci- 
ence. 


Topological Data Analysis 

In data science, one is faced with the problem of ex- 
tracting nontrivial information from datasets, which 
can ohen be represented by point clouds (i.e. hnite 
subsets) in high-dimensional vector spaces. Such 
datasets may arise from an array of sensor readings 
in an engineering testbed, from questionnaire re- 
sponses in a psychology experiment, or from popu- 
lation sizes in an ecosystem. At hrst glance, it does 
not seem like we can extract any interesting infor- 
mation from computing the homology of a discrete 
set of points. All sets of n points (with the discrete 
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▼ Pigure 3: A point cloud with a hole in it. 



topology) are topologically indistinguishable, so cal- 
culating the homology of this will count the number 
of points and give us no more information. In partic- 
ular, it cannot detect any holes. 

On the other hand, if we look at the example of a point 
cloud in Pigure 3, we intuitively see that there is a 
’hole and in particular that some nontrmal topology 
is at play. 

Without going into any concrete example, it is also 
quite clear that this hole that we see intuitively is 
a piece of interesting information about the given 
dataset that we would like to extract. Now, we need 
to formalise what we mean by this hole so that we can 
generalise it to higher dimensions and eventually get 
some information about given data sets in arbitrarily 
high dimensional spaces. There is a computationally 
feasible approach to do this, which is known as per- 
sistent homology’ 

Persistent Homology 

As we saw, it does not make sense to compute the 
homology of the given point cloud as it is. Our in- 
tuition tells us that we should draw lines between all 
the points that are close to each other, then fill out all 
the little triangles and so on to get a simplicial com- 
plex. We do not need the formal dehnition of a sim- 
plicial complex for our intuitive approach here, so 
just think about lines, triangles, tetrahedra and fur- 
ther generalised triangles glued to each other along 
their faces. It is possible to compute the homology of 
such a simplicial complex on a computer, so comput- 
ing the homology of the emerging simplicial complex 
will give us the interesting information: Our data 


set has one one-dimensional hole. This intuition is 
perfectly right, but comes with one big problem: At 
which length scale should we stop drawing lines and 
triangles? The cutoff should be large enough such 
that all the ’noise’ on a small scale closes. But if it 
is too large, we draw lines and triangles between all 
points such that the interesting hole closes and com- 
puting the homology will not give us any interesting 
information. So we have to find this cutoff dynami- 
cally. We will now sketch the idea of how to do this. 

Given a point cloud D C (which is just a fi- 
nite metric space with the induced metric of we 
imagine to draw balls of radius e > 0 around each 
point and draw a line between two points x,y G D 
Be {x) n (y) 0. We draw triangles between three 

points if the balls around any two of these points in- 
tersect each other and we proceed in the same way for 
higher dimensions. The resulting simplicial complex 
is known as the Vietoris-Rips-Complex of D at scale 
e, written YK^{D). It can formally be dehned as the 
abstract simplicial complex 

VRg(T)) \= {a (Z D \ ||x — ^ll < 2e for all x, ^ G cr} 

The following example in Pigure 4 shows the Vietoris- 
Rips-Complex of an annulus-shaped point cloud D 
for different values of e (Here, D stands for dough- 
nut and not data). 

Observe, that even though e > 0 may be varied con- 
tinuously, the Vietoris-Rips-Complex only changes 
at discrete e-values ei < • • • < e^ since there are 
only hnitely many points in D. So we arrive at a fil- 
tered simplicial complex 

YRo{D)^YR,,{D) ^ 

^YR,,{D)^---^YR,^{D), 

where we can compute the homology at each liltra- 
tion step to get 

H,{yRo{D)) ^ ^ 

^H,{YR,,{D))---^H,{YR,^{D)). 

The mathematical structure of 
TDA 

We saw above, that for each i < j we 
have a natural inclusion of simplicial complexes 
VRg. {D) ^ {D), which induces a natural map 

H^{YRe^{D)) i7>^(VRg^.(D)) by the hmctoriality 

of homology. 


48 



A Pigure 4: VRg(D) for six different, increasing e-values. 


So we have a functor from the partially ordered set 
{{0,..., n}, <) (considered as a category) to the cat- 
egory of i?-modules (in this article we only consider 
the case R = Q). This motivates the following deh- 
nition: 

Deiinition 1. A persistence module over a partially 
ordered set (P, <) with coefficients in a commutative 
ring R is afunctorfrom (P, <) considered as a cate- 
gory to the category of R-modules. 

This more abstract approach gives rise to surpris- 
ing applications of persistent homology to inner- 
mathematical disciplines such as symplectic topology 
(see e.g. [4]). In these applications, the concept of a 
barcode is crucial as well. We will discuss the intu- 
ition behind these in the following section. 

Barcodes 

It is common to depict the homology groups 
iJjt(VRg(D)) in a so-called barcode. In Pigure 6 we 
can see the barcode of the above example from hgure 


Pigure 4. Staying close to this example, we would like 
to explain how this barcode is constructed and which 
information is encoded in it. 

The barcode represents the homology classes as hor- 
izontal lines that are born at some value ei and die 
at some other value > ei. We use different colours 
for dilTerent values of /c G N. The difference e^ — ei is 
called the ’hfetime of a homology class. This barcode 
hnally gives us all of the quahtative information we 
were looking for. The homology classes with a long 
lifetime can be understood to be ’holes on a reason- 
able scale’ while the nontrmal homology classes that 
have a short lifetime (compared to the class with the 
longest lifetime) can be seen as noise in the data. The 
’important’ classes with a long hfetime can be inter- 
preted to persist over a long parameter range, which 
justihes the name ’persistent homology’ In Pigure 
6, the longest horizontal green line represents such a 
persistent homology class. This persistent homology 
class, in turn, represents the hole in our doughnut- 
shaped data set D. The other relatively long green 
lines are seen because our original data set is still quite 
far away from an ideal doughnut. If we would com- 
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pute the barcode of Pigure 3 we would see one very 
long green line and other very short green lines for 
small e-values. 

To put it in a nutshell, we can extract information 
from the barcode by drawing a vertical line at fixed 
e > 0 (the blue dashed lines in Pigure 6). Then we can 
count the intersection points with the green lines and 
get the number of holes of VRg(T)) or in other words 
/3i(VRg(i))). Counting the number of intersections 
of a blue line with the red lines gives us the number of 
connected components of VRg(T)) or in other words 
/3o(VRe(T))). The longest horizontal lines represent 
the hmportant’ topological features of the dataset. 

Now, where the formahsm is operational, we should 
argue that the above example does not demonstrate 
the full power of persistent homology. We only ob- 
served a 1-dimensional hole. In principle, we can 
compute holes of arbitrarily high dimension simul- 
taneously and compare their length scales. Think for 
example of a data set D C which approximates 
a 2-sphere with a small disk removed that we callM 
(see Pigure 5). Here, D stands for Death star not for 
Data. Here we would see no non-trivial homology 
for a large range of e (since a sphere with a ’hole cut 
out is contractible), until for large enough e we obtain 
non-trivial second homology, suggesting indeed that 
our space looks like a 2-sphere with a ’l-dimensional 
hole’ 




A Pigure 5: Roughly what M = \ D^ looks like. 
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Archimedes' Tool Kit 

Eric L Bandurski 


A skilled person of Archimedes’ time could have 
made accurate models of all 18 Platonic and 
Archimedean polyhedra by inscribing circles and 
chords on spheres. Using a few measurements from 
geometry, the tools required would be: a lathe, di- 
viders, and a drill. 

Using a knowledge of six ordinary polygons, one can 
construct all 18 Platonic and Archimedean polyhe- 
dra. One must hrst abandon the idea that mea- 
suring precise angles is necessary. A woodworker 
making molecular models soon realizes that depend- 
ing on protractors is crude. If you want precise an- 
gles, use sines, effectively chords - the old-fashioned 
way. There is a tendency of those who are per- 


haps overschooled in maths to bludgeon the prob- 
lem of constructing these polyhedra by applying ad- 
vanced mathematics. Instead, lets simphfy and let 
the beauty of rigorous geometry and the concepts of 
circle, chord, and sphere do all our work. 

Our mathematical tool kit 

The Platonic and Archimedean polyhedra use one or 
more of six regular polygons. 3, 4, or 5 regular poly- 
gons meet at each vertex. The polygons can be the 
same at each vertex (5 Platonic polyhedra) or mixed 
(13 Archimedean polyhedra). 



◄ Pigure 1 : The lengths of 
the diagonals, BE, across 
two sides of these six poly- 
gons were surely known 
to Archinnedes. He could 
easily have constructed and 
measured them using di- 
viders and straightedge. 

In Ptolemy's time (c. AD 
100--170), these chords cor- 
responding to the interior 
angles of the polygons were 
used in compiling the Al- 
magest. 
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Seeking the Simple Approach 

I knew from my reading that vertices of Platonic and 
Archimedean polyhedra are interchangeable. This 
requires that they be located on the surface of a 
sphere, which is called the circumsphere. While 
studying Wikipedia drawings I made the second criti- 
cal observation that the neighbouring vertices imme- 
diately surrounding each vertex lie on a circle. They 
must lie on a circle because they lie on the inter- 
section of two spheres: a smaller sphere whose ra- 
dius equals the edge length, S, and the circumsphere 
which is the locus of all vertices. The smaller sphere 
I will call the vertex sphere. The intersection of these 
spheres, i.e., the circle, I will call the vertex circle. 

Eureka! Our problem simplihes to precisely htting a 
set of chords of known length around a circle. This 
suggests a purely geometric method for making any 
of the eighteen Platonic and Archimedean polyhe- 
dra. We don t need a complicated apparatus for set- 
ting angles; we only need: 1) the radius of this ver- 
tex circle, 2) a chord length for scribing this circle on 
a sphere, and 3) the lengths of the diagonals for the 
regular polygons meeting at each vertex. 


let s set the radius of our vertex sphere, therefore the 
edge length of our polyhedron, S = 1. Chords of the 
four polygons must fit perfectly on the vertex circle 
of smaller radius, Tyc- 


lcosidodecahedron 



A Pigure 2: View fronn above 
vertex A of four surrounding vertices 


The icosidodecahedron 


To illustrate this concept we will construct the geom- 
etry of one of the most beautihil of the Archimedean 
polyhedra, the icosidodecahedron. To make a 
long story short (because its fully 
expanded in another paper), each of the 
18 Platonic and Archimedean polyhedra can 
be constructed precisely in ball-and- 
stick form by scribing a vertex pattern 
on accurate spheres, one for each vertex. 

The icosidodecahedron (designated 
3,5,3,5) has 30 vertices, each of them 
with four polygons (triangle, 
pentagon, triangle, pentagon) 
arranged around the vertex, as 
shown in Eigure 2. Eor convenience. 



Pigure 3: Geometry of the icosidodecahedron derived from the radius of the vertex circle A 
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1. Pinding the vertex circle radius 

The chords, BE, EF, FG, and GB (diagonals of 
polygons) must ht perfectly on the circle of radius 
Tyc- From Pythagoras we can find the radius of the 
icosidodecahedron vertex circle, r^;c, a projection of 
edge S, as ^ 0.951056. The radii of the vertex 
circles of all except two of the remaining polyhedra 
can be calculated analytically. The radii of the vertex 
circles for the snub cube and the snub dodecahedron 
can be calculated very precisely with a few iterations. 

2. Constructing a cross-section 

Once the radius of the vertex circle, has been de- 
termined, a cross-section of any of the 18 polyhedra 
can be constructed with dividers and straightedge, as 
seen in Pigure 3. Beginning at O, we lay down XB, 2 i 
diameter of the vertex circle, and Y Z, its perpendic- 
ular bisector. We locate vertex A by striking a unit 
arc from B, intersecting YZ. The same arc from A 
defines the vertex sphere and locates H. HB is an 
important chord to be used later. VW, the perpen- 
dicular bisector of edge AB, intersects YZ and de- 
fines G, the centre of the circumsphere. 

3. Using a model sphere 

Any small sphere centred at vertex A can serve as a 
perfect miniature of the vertex sphere of radius AB, 
including the angles. Pigure 5 shows a cross section of 
the model sphere. Before locating holes for the dow- 
els, chord H'B' = HB -ris used to scribe the vertex 
circle on a model sphere. 


4. Laying out chords 

Below is the model sphere viewed from G, the centre 
of the polyhedron. The chords will have lengths pro- 
portional to the smaller radius: B'E' = r, E'F' = 
Lp • r, F'G' = r, and G'B' = r. Using the dividers 
we can step off these four chords around the minia- 
ture vertex circle described earlier: first, B'E', E'F', 
and F'G', then G'B'. This should bring us back to 
B’ and complete the quadrilateral. This is a valuable 
check for precision. 


Icosidodecahedron 



A Pigure 4: Model sphere viewed fronn G 

The final step is to drill four holes of equal depth, ex- 
actly toward the centre of the sphere, at the four loca- 
tions B', E', F', and G'. Thirty of these spheres and 
sixty edges of equal length close perfectly into a rigid, 
classic icosidodecahedron. A picture is shown on the 
next page. 



A Pigure 5: Cross-section of the model sphere 
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What we used 



A Figure6 

The icosidodecahedron I constructed 


The tool kits required to make accurate ball-and-stick 
models of polyhedra turn out to be quite simple: 

Mathematical Tool Kit: chord lengths (polygon di- 
agonals), an iterative or analytical method tojind r^c 
(vertex circle radius), dindgeometry to find H'B' (arc 
length to scribe a miniature vertex circle on a sphere). 

Physical Tool Kit: lathe (spheres of equal diameter 
for vertices and dowels of equal length for edges), di- 
yiders, and a drilling apparatus (holes of equal depth, 
drilled toward the exact centre of each sphere). 

Could the ancients in Archimedes’ time have con- 
structed his 13 polyhedra, using geometry and the 
simple tools available? I believe the answer is yes. 
In fact, it appears that an approach using circles and 
chords scribed on spheres is the simplest and most 
precise way that the polyhedra can be constructed. 

Accurate spheres could have been made of a sta- 
ble wood, e.g., olive heartwood in Archimedes’ time. 
Portunately for us, precision spheres of polyethylene 
or polypropylene are used in ball valves, so theyre 
about as afFordable as wood and usually accurate to 
within a thousandth of an inch. A stringent test of 
this construction is the creation of a chiral pair of 
snub dodecahedra - my original goal: 



▲ Figure7 

A chiral pair of snub dodecahedra 
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La Baguette 
Mathĕmagigue 

Julyan H. E. Cartwright 


If you throw a needle or stick at random onto a floor 
ruled with parallel lines, such as the cracks between 
Aoorboards or tiles, from the proportion of times that 
the stick lands crossing a crack you can estimate tt; 
can we get e as well? Yes, we can. All of these aspects 
have been discussed before, but I havent seen them 
discussed in this way: that one can estimate both tt 
and e with the same BuAbns Needle experiment. 


The experiment 

This algorithm for the computation of tt by Monte 
Carlo methods is BuAbns Needle, flrst described by 
the Count Buflbn in 1733 in a report to the Acadĕmie 
des Sciences [1] and then discussed in a supplement 
to his great work of Natural History published in 1777 
[2]. For its implementation we need a floor with par- 
allel lines - the cracks between Aoorboards, or floor 
tiles if we agree to consider only the lines running in 
one direction. The lines are a distance a apart, onto 
which a pen of length l is thrown at random. The final 
position of the pen on the floor can be given in terms 
of the distance x of its midpoint from the nearest line 
and the angle (/) it makes with the lines. For the pen 
to cross one of these lines, it must be that 

X < //2 sin 0, 

and the probability of this happening is 

1 r 1 • 2/ 

p = — / Ismcpdcp = —. 

UTT Jq aTT 

This equation can be used to compute tt in the follow- 
ing way: we throw the pen at random onto the floor 


a large number of times n and count the number of 
times m it lands on a crack. Then tt is approximated 
by 

2ln 

7T ~ -. 

am 


Bread and pistols 

As the French word for stick is “baguette”, some writ- 
ers have mistakenly assumed that Buflon was advo- 
cating throwing bread around, and this error is still 
extant. In fact, Buflon himself seems neither to have 
done any needle or stick throwing, nor to have con- 
sidered his thought experiment as an experimental 
approximation to tt. Buflon was thinking of a gen- 
eral class of problems in geometric probability. 

“But if, instead of throwing into the air a round object 
such as a coin, one threw an object of another shape, 
such as a square Spanish pistole [a type of coin], or 
a needle, a stick, etc., the problem would demand a 
little more geometry, although in general it would al- 
ways be possible to give its solution by comparisons 
of space, as we are going to demonstrate”, he wrote 
[3]. 

Approximating tt 

It seems to have been Laplace who first associated 
BuAons idea with throwing a stick to obtain tt [4]; 
Laplace also solved the slightly diAerent problem us- 



ing square tiles [5, 6] (see also [7]). There are many 
more extensions: one can also consider convex bod- 
ies instead of needles [8]; or flexible bodies, like noo- 
dles [9, 10]. 

As n increases, better estimates of tt are obtained, 
albeit very slowly; since the standard error is e = 
\/p{l — p)/n, we have n ~ l/e^; around a mil- 
lion throws give three signihcant hgures. Mario Laz- 
zarini claimed in 1901 to have obtained tt to 6 deci- 
mal places by throwing a needle 3408 times [11], but 
as pointed out by Gridgeman [12] and Badger [13], 
he almost certainly cheated by choosing the best mo- 
ment to stop. We expect the estimate obtained by 
needle throwing to cross the true value of tt inhnitely 
many times as n increases, so if we truncate near one 
of the crossings we get an extremely good result, but 
we can only do this by knowing beforehand what the 
answer should be. If one doesn t cheat, convergence 
is slow and BulTons Needle is not a quick way to es- 
timate tt. 

On the other hand, it seems that nature may use the 
technique. Ants may use a Buffon-Needle-type al- 
gorithm to evaluate the size of potential nest sites 
[14]. They do this not by throwing anything, but by 
making use of the similar result that the estimated 
area of a plane is inversely proportional to the num- 
ber of intersections N between two sets of lines, of 
total lengths S and L, randomly scattered on to it: 
thus 2SL/{7vN). “Scouts [i.e., ants] using such a 
Buffons Needle algorithm will assess nest area as in- 


versely proportional to the number of intersections 
they make between a tirst set of pheromone-marked 
paths and a second set of census paths” [14]. 


Approximating e 

But we can obtain e at the same time as we get tt from 
this experiment. We get e from the result that if we 
sample random numbers uniformly from [0,1] un- 
til the sum tirst exceeds 1, the expected number of 
draws is e. Russell [15] indicates that the result comes 
from Gnedenko [16]: “an exercise in Gnedenko [...] 
asks the reader to show that, if Ui, U^, ... are i.i.d. 
uniformly on (0,1), if = Y/i=i ^ 

the minimum value of n for which Sn > l, then 
E{N) = e.” He adds, “The proof is not dihicult, and 
one approach is to solve the general problem of the 
expected number of draws for the sum to exceed x, 
withO < X < 1. You get an integro-differential equa- 
tion, the solution of which is e^.” 

So, to estimate both tt and e, throw the needle on 
floor, and at each throw, as well as recording whether 
or not it crosses a line, record its distance from the 
line nearer to you as a proportion of the distance to 
the next line — any point of the needle will do, an 
end or the middle, as long as it is consistent. Then tt 
is given by the proportion of times the needle crosses 
a line. And e is given by the average number of obser- 
vations that need to be added together to reach one. 
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Magnetic Monopoles 
in Motion 

Vit Sriprachyakul 


In any standard textbooks about quantum mechan- 
ics, we can find an introduction about the idea of 
magnetic monopoles and its constraint on values 
of charge, known as Diracs quantisation condition. 
However, an approach taken to construct a monopole 
is rather mathematical and non-intuitive. I will ap- 
proach this in yet another way, the whole idea is be- 
lieved to mainly due to Dirac. Also, I will push this 
model hirther and explore more properties of mag- 
netic monopoles. 

Throughout this article we use the following conven- 
tions: The signature is (+, -, -), the unit system 

is SI, and the Cartesian coordinate system is right- 
handed. 


Magnetic monopoles at rest 

Let me tell you what a standard method is. We start 
with a vector potential 


We then take a curl in spherical polar coordinates to 
get a magnetic field 



and by multiplying A by , we recover a point-like 

field, 

^ -O • 

47r 

Notice that A is not regular everywhere! When 
sin ^ = 0, i.e., ^ = 0 or tt, then A is singular. When 


^ = 0, this is cancelled by the numerator, however 
this is not the the case for 0 = 7v, so the negative 
2 :-axis is a genuine singularity. This is crucial since 
Gausss law says that we cannot have monopoles, so, 
we cannot have monopole field which is regular ev- 
erywhere. The singularity along the whole negative 
2 ;-axis is known as Diracs string. 

Dirac's string 

Let us explore this singularity. Consider a semi- 
inhnite solenoid along the negative ^-axis. We are 
going to compute the field it produces. 



A Pigure 1 : A senni-infinite cylindrical solenoid 

We will assume that the current is uniform and flows 
around the 2 ;-axis in a positive sense. We also assume 
that the solenoid is cylindrical and its radius is small. 
Note that if we divide our solenoid into many thin 
rings with thickness 6, then we can approximate each 
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ring as a magnetic dipole! 

The tield produced by an ideal dipole is 
^ Mo (' r)r 

^^4^ V ’ 

where m is a dipole moment; in our case, this is 
5rh = a5zAz with A a cross-sectional area and a 
a current per unit length. 

In cylindrical polar coordinates, we have 

OO / 

^ M I 3(z + zo)^ 

Jg ++ Po)^ 

1 

“((^ + ^o)2 + p§)t 

OO y 

B = / ^ I 3(z + zo)po 

Jo ^^[((^ + ^o)^ + Po) 

where (zo, Po? ^o) is the point at which we evaluate 
the field and 2 ) is the position of our thin ring, mea- 
suring downwards. 

This can be easily integrated. If you do not want 
a hint, skip the next sentence: Use the substitution 
z + zo = po tan h, or by parts. 




This is exactly a monopole field! We then identify 
magnetic charge, qmy with aA, which is a flux inside 
the solenoid divided by po- We can show that mag- 
netic charges are Lorentz invariant, just like electric 
charges. This is one of the topic we are exploring be- 
low. 


...in motion 

Let us take magnetic charges to be the solenoid we 
discussed above. We see that the charge is cor- 
responding to the flux through the interior of the 
solenoid. Let S be an inertial frame that the solenoid 
is at rest and 5" be an inertial which sees S fly off with 
a constant velocity u in their mutual x-direction. 


s 


S' 







▲ Pigure 2: The two trames of reference 


We then end up with 


B,= 

Bp = 


+0 

47r 

47r 


aAzo 

(^o + Po) ^ 

crApo 

(^o + Po) 


For the magnetic charge to be invariant, we need to 
prove that the flux through the solenoid is the same. 
We can do this by looking at how the magnetic and 
electric transform under the change of frames. Note 
that the solenoid can lie in any direction, with its end 
fixed at the origin as shown. 
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We have 

E' = —^uU X B. 

We do not need the second equation for now, but it 
is relevant for later. We see that if we choose the sur- 
face for evaluating the flux to lie in the mutual Y — Z- 
plane, this surface does not have a length contraction 
effect, since it is perpendicular to the relative veloc- 
ity. Since the field in the x direction is also the same, 
we conclude that the flux is the same in both frames 
and thus is Lorentz invariant as claimed. 

There are some subtleties to be explained, such as 
the question “Is the flux through our clever surface 
the same as through the surface normal to solenoid 
axis?”. This is easy: since the magnetic field inside the 
solenoid lies in the same direction as the solenoids 
symmetry axis, the flux through the two surfaces are 
the same. 

Secondly, what if the solenoid lies in the Y—Z -plane? 
Now we cannot use our clever surface, but this case is 
actually simple! We choose our surface to be the one 
normal to the symmetry axis. This surface lies in the 
X — Z-plane and thus suffers a length contraction ef- 
fect in the x-direction, so the total cross sectional area 


is cr' = But the Aeld in the ^-direction also trans- 

lu ^ 

forms By = ^uBy. These two factors of ^u cancel 
and we still have the invariant property of magnetic 
charge! 

Now in the velocity limit \u\ c, we can use the 
transformation for E to see that, in the frame that 
sees the charge move. 


This is the Biot-Savart Law for a magnetic charge. 

Interaction 

For our hypothetical magnetic charges to have a 
chance to be real physical charges, the above Lorentz 
invariant property is expected. Otherwise, one can 
send the monopole into motion and its charge trans- 
forms, so we know it is not hmdamental. It could 
even break Diracs quantisation as well, since the 
Lorentz transformation is continuous and if the 
charge did transform, it should do so in a continu- 
ous fashion. Dirac s condition however is discrete, so 
it is bound to break the quantisation condition. 
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Nonetheless, this is not the only thing we can assume: 
we can expect to derive the interaction with electro- 
magnetic helds as well, which is our main interest in 
this section. 

Suppose we have two magnetic charges, both at rest, 
what is the force between them? We take one of the 
charge to lie in negative 2 ;-direction and this subject 
to a monopole held from the other. We again slice 
our solenoid into thin rings. 

The force acting on a dipole is given by 

F = V{m ' B). 

Notice that the gradient is differentiating w.r.t. the 
position of the dipole and this is inconvenient since 
we need to sum all these tiny forces to get a resul- 
tant force and the derivative is differentiating differ- 
ent points! Luckily, in this case of a monopole held, 
one can show that V = where R is the po- 

sition vector of the monopole producing the B held. 
The total force is then 

F = — B ' dm 

So we only need to compute the following integral: 

r {zFzo)aAdz 

J 47r + + 

dTTR 
So 

1 2 

m _ I^OQmQm £> 

This is Coulombs law for magnetic charge. No- 
tice that same-signed charges repulse and opposite- 
signed charges attract. Note also that the superscripts 
here are indicators for the charges, not exponents. By 
arguing that general magnetic helds are a superposi- 
tion of monopole helds, one can deduce that, for a 
general magnetic held, 

F = qmB. 


j B • dm = — 


The expected Lorentz force law is then 



The above derivation is approximated, so how can we 
trust this result? Well, at least it gives us a hint of what 
to write down! Since we need our physics to be con- 
sistent with special relativity, we should have 




dp^ 

dr ’ 


where (i.e. is the Hodge dual of 

is the 4-velocity of the magnetic monopole, 
and = mu^ is the 4-momentum. 

If you expand the above equation, you should recover 
the Lorentz force law stated previously. 


Dirac's quantisation condition 

There is no detailed calculation in this section, only 
an outline of how the condition can be derived is pre- 
sented. 

We can take our solenoid model of monopoles, and 
then try to observe its presence by using the quantum 
double-slit interference experiment. We can calcu- 
late the variation in probability amplitude from the 
Peynman path integral and the Lagrangian for a par- 
ticle in an EM held. 

We then say, for the monopole to exist, we need the 
solenoid to be not observable. This is because if we 
can observe the solenoid, then we know that the held 
is not everywhere point-like and thus not a genuine 
point charge. This gives a constraint on quantum in- 
terference amplitude and so gives Diracs quantisa- 
tion condition. 


Next, we are interested in how it interacts with an 
electric held. Lets assume for simplicity that we are 
in a low speed limit. Suppose that in our frame, the 
monopole is moving with velocity u in a background 
EM held E^B. By transforming these helds to an in- 
stantaneous rest frame of the monopole, we get 
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Makers of Patterns: 
From Escher to Coxeter 

Alexis Marchand 


Dutch artist M. C. Escher was fascinated with the idea 
of representing inhnity within a hnite piece of art. 
Among many sources of mathematical inspiration, a 
correspondence with the geometer Donald Coxeter 
led him to his Circle Limit series, an example of which 
is given in Eigure 1. Unsurprisingly, this series of 
drawings hides deep mathematical ideas connected 
with the work of Coxeter, and we are going to see how 
these ideas can lead to modern concepts in geometric 
group theory. 

Examining Circle Limit 7, we see that it is constructed 
by starting with a simple building block, a (hyper- 
bolic) polygon with some decoration (shown in Eig- 
ure 2), and then rehecting this initial polygon along 
its edges, and iterating the process. 



A Pigure 1 : Cirde Limit i, M. C. Escher. 



A Pigure 2: TheEundannental domain 
of Eschehs Cirde Limit i 


In order to formalise this construction, we shall work 
in a space X that can be either the euclidean plane E^, 
the sphere or the hyperbolic plane H^, the latter be- 
ing the setting of Escher s Circle Limits. A geodesic in 
X is the shortest path between two points: euclidean 
geodesics are straight lines and spherical geodesics 
are great circles. Any geodesic can be extended in 
both directions to either an inhnite geodesic (in the 
euclidean and hyperbolic cases) or a closed geodesic 
(in the spherical case); there is then an operation of 
rehection with respect to that geodesic, which dehnes 
an isometry of X. The hyperbolic plane can be repre- 
sented, as in Escher s work, by an open disk in which 
geodesics are arcs of circles (or straight lines) orthog- 
onal to the boundary; this representation is called 
the Poincarĕ disk model. In this model, distances are 
more and more distorted as one gets closer to the 
boundary. One can think of it as a way to visualise 
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a geometric universe which does not obey euclidean 
rules, in the same way that a standard Mercator pro- 
jection map of the world will alter distances and mag- 
nify the regions that are closer to the poles. 

We consider a compact convex polygon P in 
X, in other words a set of vertices together 
with geodesic segments between them; this is 
our tilings basic building block. We denote by 
xo, xi,..., Xn-i,Xn = xo the successive vertices 
of P, and we assume that the interior angle at 
can be written as n/pk for some integer > % 
as in Pigure 3. For instance, the hyperbolic poly- 
gon of Pigure 2 has four vertices, with interior angles 
tt/S, 7r/2, tt/S, 7r/2. 


X3 



A Pigure 3: A polygon satistying the hypotheses 
ofthe Poincarĕ Theorem 


Let : X —)► X be the rehection of X with re- 
spect to the geodesic segment [x/g_i, x/g]. The tiling 
we are interested in can be constructed by drawing 
all the polygons obtained by applying sequences of 
rehections among the cr/gS to the original polygon P, 
or in other words, by making the group L generated 
by CTi,..., cr^ act on P; this group L is therefore our 
hmdamental object. 

The Poincarĕ Theorem 

Because L determines the tiling, we wish to under- 
stand its algebraic structure. We would like in partic- 
ular to hnd a presentation of T : this means hnding a 
small set of algebraic relations between the generators 
that determine all other relations that exist in T. We 
start by writing the most obvious relations we have: 
since is a rehection, it satishes cr^ = 1. Moreover, 
one can check that the composite o cr/g_^i is a rota- 
tion around Xk of angle ln jpk\ therefore it has order 


p/c, which implies that (cTj^crj^+i)^^ = 1. Nowitturns 
out that these relations entirely determine L, as stated 
by the following theorem. 

Theorem 1 (Poincarĕ). Thegroup L = (cri,..., a^) 
has the/ollowing presentation: 

r ^ (cTi,..., cr„ I Vfc, {akCTk+if‘ = 1, al = l). 

( 1 ) 

Moreover, T acts properly on X and P is afundamen- 
tal domain for this action (i.e. the T-translates of P 
cover X). 

Saying that T has a presentation given by (1) means 
that r is, in some sense, the simplest group gen- 
erated by cri,..., cr^ and satisfying the given re- 
lations. More precisely, T is isomorphic to the 
free group on {cri,..., cr^} quotiented by the nor- 
mal subgroup generated by {cr^, l <k <n} U 
5 1 < ^ < 'ki}. The fact that T acts 
properly on X means that the action is, in some sense, 
discrete: in the resulting tiling, there cannot be an in- 
hnite number of copies of P in a small region. 

The proof of the theorem involves reconstructing the 
tiling as a cell complex, i.e. a space obtained by glue- 
ing various cells. This cell complex K is dehned by 
taking one copy P^ of P for each element 7 of T and 
by glueing the k-th edge of P^ with the k-th edge of 
Pjak all 1 < /c < n. One has to show that K is 
homeomorphic to X and that the action of T on K 
corresponds to the action on X, and this allows one 
to actually understand the structure of T. For more 
details, the reader is referred to [4]. 

Before going any hirther, let us examine a few spe- 
cial cases of the Poincarĕ Theorem. If we look back at 
the polygon of Pigure 3, we see that it is impossible in 
the euclidean world: the sum of angles of an n-sided 
polygon is (n — 2 ) 77 . We should therefore have 



But since > 2 for all k, the right-hand side of the 
above equation is at most n/2, which is only possible 
if n = 3 or 4. Reasoning along these lines, we see 
that the only euclidean polygons to which the the- 
orem applies are rectangles, the equilateral triangle, 
and right-angled triangles with angles 7r/2,7r/3,7r/6 
or7r/2,7r/4,7r/4. Pigure 4 shows such a tiling. 



A Pigure 4: Tiling of by a triangle with angles 7r/2,7r/4,7r/4 


In the spherical world, the sum of angles of a trian- 
gle is always greater than tt; as a consequence, the 
only polygons satisfying the hypotheses of the the- 
orem are triangles, but there are inhnitely many of 
them, for instance with angles 7r/2,7r/2,7r/p for all 
p > 2. In the hyperbolic world on the other hand, 
because the sum of angles of a triangle is less than tt, 
there is a lot more freedom. For instance, hyperbolic 
regular right-angled n-gons exist for all n > 5. Fig- 
ure 5 shows both a hyperbolic and a spherical tiling. 



A Pigure 5: Tiling of by a right-angled hexagon and 
of by a triangle with angles 7r/2,7r/2,7r/5 


Coxeter groups and 
Cayley graphs 

The group T with the presentation given by the 
Poincarĕ Theorem is part of a class of groups intro- 
duced by Donald Coxeter in a 1934 article [1] in 
which he undertook a study of “discrete groups gen- 
erated by reHections” in euclidean and spherical ge- 
ometry; these groups would later become known as 
Coxeter groups: a group TP is a Coxeter group if it 
has a presentation of the form 

w = {S S, = {st)^‘* = 1) , ( 2 ) 

where 5 is a hnite set and (mst)g t<£S ® sym- 
metric matrix indexed by S, with coethcients in 
{2,3,...}U{oo}. 

The simplest example of a Coxeter group is 
the one with only two generators: W = 
(s,t I = (st)^ = l). In order to visualise 

this group, let us introduce the concept of Cayley 
graph: given a group G with a generating set S, its 
Cayley graph Cay(TC, S) is the graph obtained by 
drawing one vertex for each element of G, and by 
drawing an edge labelled by s between g and gs for 
all g C G. Hence, the local model for a Cayley graph 
is the foIIowing: 

s 

9 • •gs 


Coming back to our previous example of 
W = {s,t \ s^^ = = (st)'^ = l), we can enumer- 

ate elements oiW: starting with the trivial word 1, 
we obtain s, st, sts, etc. These words are all dif- 
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ferent until we reach = 1. We could also 

have started with t, ts, etc., but we would have ob- 
tained the same words because t = {st)^~^ s, ts = 
{st)'^^^, etc. Therefore, if m is hnite, the Cay- 
ley graph Cay {W,{s,t}) should look like a cycle of 
length 2m, with edges labelled by s and t alterna- 
tively. 



A Pigure 6: The Cayley graph of 
Dq = (^s,t \ s'^ = t^^ = {st)^ = l) 


If m = oo, we can multiply indehnitely by s or t, so 
the Cayley graph of W is an inhnite line. 


tst ts t t 1 S S St sts stst 

k Pigure 7: The Cayley graph of 
Doo = {s,t \ s^^ = t^^ = l) 


It turns out that having determined the Cayley graphs 
of Coxeter groups with two generators tells us how 
to draw the Cayley graph of any Coxeter group: given 


A Pigure 8: The Cayley graph of 
(r, s, f I = (rs)^ = {rt)^ = (st)^) as 

the dual of the corresponding Euclidean tiling 


W as in (2), if one looks at edges labelled by two gen- 
erators s,t ^ S only in Cay(W, S), then one will see 
cycles of length 2mst (or lines if mst = oo). 

If we look back at the cell complex K that we used 
to reconstruct the tiling of the Poincarĕ Theorem, 
we will see that we have done something very simi- 
lar to the construction of Cayley graphs: we started 
with a collection of polygons (instead of vertices for 
the Cayley graph) indexed by elements of the group, 
and we glued edges of the polygons (instead of link- 
ing vertices by an edge) when the two correspond- 
ing elements of the group differ by a generator. This 
suggests the following construction illustrated in Fig- 
ure 8: in the tiling, add a vertex inside each copy of 
the polygon. If two copies of the polygon share an 
edge, then draw a new edge between the correspond- 
ing vertices across the old edge. The resulting graph 
(consisting of new vertices and edges only) is called 
the dual graph of the tiling, and is exactly the Cayley 
graph of the corresponding Coxeter group. 

Quasi-isometries 

In the context of the Poincarĕ Theorem, the construc- 
tion of the Cayley graph as the dual of the tiling seems 
to point to a connection between the geometry of 
X and that of Cay(r, S). One way to put it is that, 
forgetting the local details and looking at them from 
afar, those two objects have very similar geometric 
structures. We can make this intuition rigorous by 
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introducing the concept of quasi-isometry. Given two 
metric spaces X and Y, a map f : X ^ Y is said to 
be a quasi-isometric embedding if there exist £ > 0 
and A > 1 such that 

jdx {x,x') - e <dY {f{x),f{x')) 

< Xdx {x, x') + 

for all x,x' G X. We say that / is a quasi-isometry 
(and that X and Y are quasi-isometric) if in addition 
/ has a quasi-inverse: a map g \ Y ^ X such that 

sup dx {x,go f{x)) < oo 

xex 1 

and 

sup dy {y, f O g{y)) < oo. 
yeY 

Quasi-isometry dehnes an equivalence relation be- 
tween metric spaces and captures the idea of “looking 
the same from afar”. 

Given a group G with a generating set Sy the Cayley 
graph Cay(G, S) can be made into a metric space by 
endowing it with the graph distance: the distance be- 
tween any two points in the graph is the length of the 
shortest path between them. If we restrict this dis- 
tance to G, we obtain the word metric ds onG, which 
can be dehned equivalently by saying that ds^g, h) is 
the length of the shortest word on S equal to h~^g. 
In fact, restricting our attention to elements of G or 
considering the whole Cayley graph does not matter: 
(C, ds) is quasi-isometric to Cay(C, S). Now one of 
the reasons for the importance of quasi-isometry is 
the following fact. 

Proposition l» IfG isagroup and S, S' are twojinite 
generating sets for G, then (G^ds) is quasi-isometric 
to {G,ds'). 

Therefore, even though a hnitely generated group 
may have many different Cayley graphs and word 
metrics, the proposition implies that they are all 
quasi-isometric. This allows one to talk of the ge- 
ometric structure of a (hnitely generated) group, 
which is well-dehned up to quasi-isometry and can 
be visuahsed as the structure of any Cayley graph of 
the group. 

The geometry of groups 

The idea of dehning the geometric structure of hnitely 
generated groups and accepting to look at metric 
spaces up to quasi-isometry was hrst introduced by 
Mikhail Gromov in the eighties; this is the hrst step 


into the realm oigeometric group theory, a very active 
held of contemporary mathematical research. One 
hmdamental result of the held is the following: 

Lemma 1 (Svarc-Milnor). Let G be a group acting 
properly, cocompactly and by isometries on a proper 
geodesic space X. Then: 

1. G is jinitely generated. 

2. G is quasi-isometric to X. 

It follows from the Svarc-Milnor Lemma that the 
Coxeter group T appearing in the Poincarĕ Theorem 
is quasi-isometric to the space X. Amazingly, we can 
use this fact to deduce many algebraic properties of 
r: for instance, if X is the hyperbolic plane, then T 
is a so-called hyperbolic group, which implies for ex- 
ample that it does not contain a subgroup isomorphic 
to (the intuition behind this is that is flat, so it 
cannot be embedded into the hyperbolic world). 

This short article will hopehilly have pointed to the 
idea that viewing groups as intrinsically geometric 
objects and trying to understand them as the sym- 
metries of some spaces is both a natural and beauti- 
ful approach. A good place to learn more about var- 
ious aspects of this general idea is [4]; a more linear 
textbook-style approach is given by [5]. For more on 
Coxeter groups, we recommend [3] for their combi- 
natorial theory and [2] for their geometry. 
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Simpli^ying 
Differential Eguations 
with Group Work 

Christopher J. Lang 


As any mathematician will testihy, hnding all solu- 
tions to a diiferential equation is a daunting task. One 
made much easier when working in groups. No, not 
a collection of people — though that certainly would 
help — rather, the algebraic structures. 

Typically, solutions to dilferential equations depend 
on several independent variables, like the initial con- 
ditions. As such, the solution space is often mas- 
sive. However, we can create equivalence classes, via 
a group action, reducing this space to a more manage- 
able size. By choosing specific representatives from 
these classes, we can solve the difterential equations 
more easily. Using our group action on our represen- 
tatives, we then recover all solutions to the difteren- 
tial equations. Moreover, multiple group actions can 
be used in conjunction to reduce our space and sim- 
phfy our task of solving these difterential equations 
even hirther. 

While this technique is general, we apply it to the 
Nahm equations: a system of first-order, coupled, 
nonlinear difterential equations. By focusing on this 
example, we can see the power of this technique first- 
hand. This article presents and elaborates on two of 
Andrew Dancers papers [ 1 , 2 ]. 

Suppose we have N, a natural number, and 
(a, 6 ), an open interval of the real line. Let 


To,Ti,T 2 ,T 3 : (a, 6) ^ ti(A^) be analytic, where 
u(A^) is the set of 7V x N, skew-Hermitian, complex 
matrices. The Nahm equations are given by 

T^ + [To,Tj] = [n,Ti], ( 1 ) 

for all cyclic permutations of (1,2,3), with 

the commutator [A, B] := AB — BA. Let r be the 
set of quadruples (To, Ti, T^, T 3 ) = T of matrix val- 
ued hmctions satisfying the following: 

1 . To, Ti, T 2 , and T 3 must be u(A^)-valued and 
analytic; 

2 . To is analytic on [a, h] and Ti, T 2 , and T 3 are 
analytic on (a, 6) with simple poles at t = a,b; 

3. T satishes the Nahm equations; 

4. the residues (from complex analysis) of Ti, T 2 , 
and T 3 at t = a,b form an irreducible repre- 
sentation of 5 u( 2 ). 

Elements of r are called Nahm data. 

The Nahm equations are obtained by imposing trans- 
lation invariance in three dimensions on the four- 
dimensional, self-dual, Yang-Mills equations. Via the 
ADHMN procedure, Nahm data is used to create 
magnetic monopoles. 
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The Power of Group Work 

A system of coupled, nonlinear difFerential equa- 
tions, the Nahm equations are extremely dithcult to 
solve. However, we can dehne group actions that al- 
low us to greatly simplify the form our solutions can 
take. This is done in the next section. 

Theorem 1. Every element (To,Ti,T^,T3) G r is 
equivalent, via the actions of the gauge group andM^, 
to Nahm data oftheform 

To = 0; Ti,T 2 ,T 3 5u(iV)-valued, (2) 

where 5u{N) is the set oftraceless matrices in u{N). 

Before using group actions, we were tasked with hnd- 
ing real, analytic hmctions. The above theorem 
ensures that we only have to hnd 3N^ — 3 such func- 
tions. In the N = 2 case, this means we have to find 
nine of the original 16 hmctions. By using an addi- 
tional action of SU(2) on the N = 2 Nahm data, we 
can reduce this to three functions, though we do not 
cover thathere [2]. 

The above theorem also has ramihcations for the 
Nahm equations themselves. Indeed, given Nahm 
data of the form in the theorem, the Nahm equa- 
tions become Tj = [T/j, T/], for all cyclic permuta- 
tions (j, /c, /) of (1,2,3). Not only have we removed 
the need to solve for N^^ + 3 hmctions, we have also 
simplihed our system of equations through our use 
of group actions. 

For T, T G T, we say T ~ T if there are elements 
of and the gauge group which act on T to give T. 
This is an equivalence relation. Our theorem tells us 
that for every T G r, there is a T equivalent to T 
of the desired form. If we solve the Nahm equations 
for T, then we can use the group actions to recover 
our initial Nahm data. In fact, we can determine all 
Nahm data in that class through application of the 
group actions on T. Hence, solutions of the form in 
the theorem generate all Nahm data, so we need only 
focus on them when solving the Nahm equations. 

Lights, Camera, Actions 

For A G and T G r, the action of on r is given 
by 

A.T := (To, Ti — i\il, — iA^/, T3 — iA^/), ( 3 ) 

where / is the N x N identity matrix. We can con- 
firm this is a group action by directly checking the 


dehnition. Note that the changes to the j = 1,2,3 
components do not afiect the residues. 

Before we use our new action, we first prove a helphil 
feature of the Nahm data. 

Lemma 1. For T G r, Tr(Tj) is constant, for j = 
1,2,3. 

Proof As the trace operator is linear, (Tr(Tj))^ = 
Tr(Tj). Since T G r, it satishes the Nahm equations, 
so Tj = [Tk,Ti] — [To, Tj]. Hence, we have to take 
the trace of commutators, which is zero. □ 

We use this lemma together with our action to sim- 
plify the j = 1,2,3 components of the Nahm data. 

Lemma 2. Up to the action o/R^, Tj is 5u{N)-valued, 
for j = 1, 2,3. That is, there exists A G R^ such that 
forj = 1,2,3, {\-T)j is5u{N)-valued. 

Proof. Suppose T e r. Consider the hmction \j := 
^Tr(Tj), which is R-valued, as Tj is u(iV)-valued. 
By Lemma 1, \j is constant, so A G R^. Hence, 
A.T G T. As the trace of / is N, Tr{{\.T)j) = 0, 
for j = 1,2,3. □ 

Let ^ be the group of analytic hmctions g: [a,b] 
//(//), where//(//) is theset ofNxN, unitary, com- 
plex matrices. This is a group under pointwise matrix 
multiplication, called the gauge group. For g ^ G 
and T ^ r, the action of ^ on r is given by 

g.T := 

(4) 

where we multiply pointwise. Again, we can conhrm 
this is a group action by directly checking the defini- 
tion. Note that as g is analytic on [a, /], we still have 
simple roots for the j = 1,2,3 components. This ac- 
tion conjugates the residues at the endpoints, giving 
a new representation, equivalent to the original. 

We use our gauge group action to simplify the zeroth 
component of the Nahm data. 

Lemma 3. Up to the action ofG, Tq = 0. 

Proof. Suppose T G r. We want an analytic,//(A/")- 
valued hmction g satisfying g' = gT^. As To is ana- 
lytic and u(A^)-valued, such a hmction exists, for in- 
stance, g{t) = exp ^ To{F)dt^. By the dehnition 
of the gauge group, g ^ G- We know g.T G r and, 
by construction, {g.T)o = 0. □ 
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Using our actions, we now prove the theorem. That 
is, we show that there are A G and g ^ Q such 
that X.{g.T) G r possesses the desired properties: 
(A.(^.T))o = 0 and {X.{g.T))j is su^A/^^-yalued for 
J = 1,2,3. 

Proof of Theorem. Suppose T G r. By Lemma 3, 
there is a ^ G ^ such that g.T G r satishes {g.T)o = 
0. Purthermore, by Lemma 2, there is a A G such 
that (A.(p.T))j is 5u(iV)-valued, for j = 1,2,3. Fi- 
nally, the action of does not atFect the zeroth com- 
ponent of p.T, so (A.(p.T))o = 0. □ 
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